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The magnetic susceptibilities \ versus temperature T of powders and single crystals of the 
ambient-pressure (AP) and high-pressure (HP) phases of (VO)2P2C>7 are analyzed using an ac- 
curate theoretical prediction of x(T, Ji, J2) for the spin-1/2 antiferromagnetic alternating-exchange 
(Ji, J2) Heisenberg chain. The results are consistent with recent models with two distinct types of 
alternating-exchange chains in the AP phase and a single type in the HP phase. The spin gap for 
each type of chain is derived from the respective set of two fitted alternating exchange constants and 
the one-magnon dispersion relation for each of the two types of chains in the AP phase is predicted. 
The influences of interchain coupling on the derived intrachain exchange constants, spin gaps, and 
dispersion relations are estimated using a mean-field approximation for the interchain coupling. 
The accuracies of the spin gaps obtained using fits to the low-T x(T) data by theoretical low-T 
approximations are determined. The results of these studies are compared with previously reported 
estimates of the exchange couplings and spin gaps in the AP and HP phases and with the magnon 
dispersion relations in the AP phase measured previously using inelastic neutron scattering. 

PACS numbers: 75.50.Ee, 75.40.Cx, 75.30.Et, 75.10.Jm 



I. INTRODUCTION 



A resurgence of research on the magnetic properties 
of low-dimensional quantum spin systems has occurred 
over the last decade. This work was mainly initially 
undertaken to understand the relationships between the 
magnetic properties of layered cuprates containing Cu +2 
spin-1/2 square lattice layers and the high superconduct- 
ing transition temperatures of the doped materials.!!] This 
goal has also spawned much research on related one- and 
two-dimensional (ID and 2D) spin systems. Indeed, the 
new subfield ofjSgjin ladder physics was created as a result 
of these efforts. B'B The basic n-leg spin ladder consists of a 
planar arrangement of n parallel vertical spin chains (the 
ladder "legs") with (nonfrustrating) horizontal nearest 
neighbor couplings between adjacent chains, i.e., across 
the ladder "rungs". 

Self-doped two-leg spin ladders are realized in the com- 
pound NaA^Os, in which the V atoms are crystallograph- 
ically equivalent and the oxidation state of the V atoms 
+4.5, resulting formally in a mixed-valent d 05 system 
However, the material is a semiconductor rather than a 
metal. Theoretical studies have indicated that the rea- 
son for this is that the one d electron per two V atoms 
is localized on the respective V-O-V rung of the ladder 
due to the on-site Coulomb repulsion on each rung. This 
in turn led to the hypothesis that each rung acts as a 
spin-1/2 site, in which case the compound is expected to 




FIG. 1. Sketch of the basic structures of AP-(VO) 2 P 2 7 
and HP-(VO)2P207 in the a-c plane, showing the ex- 
change constants J 1 and J2 along the spin-1/2 V +4 alter- 
nating-exchange chains and J a in a perpendicular direction 
along the legs of the structural two-leg ladders (adapted 
from Ref. [ll]). The lattice parameters of AP-(VO)2P207 are 
a « 7.74, b « 9.59, and c » 16.59 A, and of HP-(VO) 2 P 2 7 
are a » 7.58, b ~ 9.55, and c ~ 8.36 A. 

behave magnetically like a spin S — 1/2 antiferromag- 
netic (AF) Heisenberg chain. Additionally, the obser- 
vation of a spin dimerization transition below 34 K, ac- 
companied by a lattice distortion, is consistent with this 
scenario. These and other aspects of the thermal and 
magnetic behaviors of NaV20s were recently examined in 
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detail in a combined theoretical and experimental study.El 
As part of this study, an accurate function was generated 
for the theoretical magnetic susceptibility x versus tem- 
perature T of the spin 5 = 1/2 AF alternating-exchange 
( Ji and Ji with 3 2 < J\) Heisenberg chain over the entire 
range < a < 1 of the alternation parameter a = J%/ J%. 

The availability of this high-accuracy theoretical 
x(k-p,T / Jx, a) function now allows one to accurately and 
precisely test the consistency of proposals for the oc- 
currence of alternating-exchange chains in specific com- 
pounds by comparing the observed x(T) with that ex- 
pected theoretically. In particular, in the present work 
we use this method to test the consistency of recently 



proposed alternating-exhange chain models for the 
ambient- and high-pressure phases of vanadyl pyrophos- 
phate, (VO)2P2C>7, and we obtain the exchange con- 
stants and spin gaps in the respective chains. The in- 
fluence of interchain coupling on the derived intrachain 
exchange constants and spin gaps is investigated using 
a mean-field approximation for the interchain coupling. 
The accuracies of the spin gaps obtained using fits to the 
low-T x(T) data by theoretical low-T approximations are 
determined. The results of these studies are compared 
with previously reported estimates of the exchange cou- 
plings and spin gaps in the two phases and with the dis- 
persion relations measured for the ambient-pressure 



TABLE I. Spectroscopic splitting factor (g- factor) g, exchange constants Ji and J2 and spin gap A determined by the method 
listed for the ambient-pressure form of (VO)2P2C>7 [AP-(VO)2P20t], for the high pressure form HP-(VO)2P207 and for the V 
dimer compound (VO)HP04 ■ ^H20. Quantities marked with a daggar (+) are derived using the S = 1/2 alternating-exchange 
chain model whereas those marked by an asterisk (*) are derived using the S — 1/2 two-leg ladder model. Quantities with 
no mark make no assumption about the specific model and/or are determined or assumed independently of the fits used to 
obtain the other listed quantities. All measurements were carried out on polycrystalline samples except as otherwise noted. 
The listed alternating exchange constants from Ref. ^ were obtained within a 2D coupled alternating-exchange chain model 
with interchain couplings given in the text and were determined from a fit to one-magnon inelastic neutron scattering data 
(Ref. [ll|) with the spin gap listed and from high-T x(T) data (Ref. |l^). The error bars on the parameters listed from the 
present work include the estimated influences of interchain interactions. 
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phase by inelastic neutron scattering. 

The history of the study of the magnetic properties of 
(VO)2P2C>7 is interesting and extensive. In the following 
two sections we give brief overviews of the previous work 
on the ambient-pressure (AP) and high-pressure (HP) 
phases of this compound, respectively, and then present 
the plan for the remainder of the paper. 

A. AP-(VO) 2 P 2 7 

The V +4 d 1 ambient-pressure phase AP-(VO)2P2C>7, 
sometimes abbreviated in the recent literature as 
"VOPO", is an industrial catalyst for the selective oxida- 
tion of rt-butane to maleic anhydride B u This compound 
was found to have an orthorhombic crystal structureo'EI 
(containing four inequivalent types of V atoms) which 
can be viewed cwstallographically as containing S — 1/2 
two-leg ladders,E3 where the rungs of the ladder lie along 
the c-axis and the legs, are oriented along the a-axis, 
as sketched in Fig. [ljl!!l A single crystal x-ray diffrac- 
tion structural study claimed that the previous struc- 
tural studies were incorrect and that the structure is 
monoclinic with eight inequivalent V atoms in the unit 
cell, although the overall structural features and the unit 
cell dimensions were found to be very similar to those 
of the previously proposed orthorhombic structure!^ 
However, a recent study of a polycrystalline sample us- 
ing both x-ray and neutron diffraction Rietveld refine- 
ments and transmission electron microscopy confirmed 
the orthorhombic structure and ruled out the mono- 
clinic structure ;t3 it was suggested that whether the or- 
thorhombic or monoclinic structure occurs in a particu- 
lar sample may depend on the exact composition and the 
details of sample synthesis.ll3 There have been two con- 
ventions used in the literature for designating the b and c 
axes in which these two axes are mutually interchan ged ; 
we will adhere to the convention in Fig. [y, as in Ref. [T_|, 
for which the approximate lattice parameters are listed 
in the figure caption. 

The magnetic susceptibility versus temperature x(T) 
of AP-(VO)2P207 was found to exhibit an energy gap 
("spin gap") for magnetic excitations]!!] The x(T) was 
initially fitted to high precision by the prediction for 
the S = 1/2 AF alternating-exchange Heisenberg chain, 
with the exchange constants J\ and Ji .and alternation 
parameter a = J2/J1 given in—Table |£j A fit by the 
spin ladder model as suggestedE3 from the crystal struc- 
ture (see Figs. [j] and ||) was not possible at that time 
(1987) due to lack of theoretical predictions for x(T) 
of this model. It was also speculated that even though 
the crystallographic features suggest a spin ladder model, 
the actual magnetic interactions might turn out to ccfc 
respond to those of alternating-exchange spin chains S3 
When x(T) calculations for the spin ladder model were 
eventually donejiil itr-was found that the same experi- 
mental x(^) data setl!9 could be fitted by the spin lad- 



der model to the same high precision as far the very 
different alternating-exchange chain modelo The exis- 
tence of a spin gap was subsequently confirmed and its 
value estimated from NMR measurements of the 31 P 
Knight shift 31 if£Xl-and nuclear spin-lattice relaxation 
rate 31 (l/Ti)(T)|LaqUrom X (T) (Ref. 0) and inelas- 
tic neutron scatteringcS measurements on polycrystalline 
samples at low temperatures, i-apd most recently from 
elastic constant measurementsliia on a single crystal in 
pulsed magnetic fields up to 50 T at T = 1.6 K and from 
Raman scattering intensity, measurements on single crys- 
tals at low temperaturesj23 as listed in Table |. 

The above inelastic neutron scattering measurements 
on polycrystalline AP-(VO)2P207 reportedly confirmed 
the spin ladder model and rejected the alternating- 
exchange chain model by a comparison of the observed 
spin gap [43(2) K] with the values 45. 7K and 57K pre- 
dicted for the two respective models from the sets of 
exchange constants determined from respective fits to 
the x(T) data.tS However, subsequent inelastic neutron 
scattering results on polycrystalline samplestil and espe- 
cially on a,-eollection of about 200 oriented small sin- 
gle crystalsliil proved that AP-(VO)2P2C>7 is not a spin- 
ladder compound. The strongest dispersion of the one- 
magnon spectra of the single crystals was found to be 
along the c axis, i.e., in the direction of the structural 
ladder rungs, and the coupling in the direction of_the lad- 
der legs was found to be weakly ferromagnetic]!^ Thus, 
perhaps surprisingly, the superexchange coupling path 
V-O-P-O-V along the c axis, coupling the structural two- 
leg ladders as shown in Fig. g, is much stronger than the 
shorter V-O-V coupling along the ladder legs parallel to 
the a-axis. These results were interpreted in terms of an 
alternating-exchange chain model with the chains run- 
ning along the c axis, with weak coupling between the 
chains, thus confirming the above speculation in Ref. [LQ- 
Subsequent x(T) data for powderE 3 ] and single crystalEHI 
samples have been interpreted in terms of the alternating- 
exchange chain model, with exchange parameters listed 
in Table | 

The spin gaps of 43(2) K (Ref. [l|) and 40.4(4) K 
(Ref. |2l|) found from inelastic neutron scattering mea- 
surements on powder samples are both significantly 
larger than the one-magnon spin gap of 36.2(3) K de- 
termined frorrL-the neutron scattering measurements on 
single crystals.!!!! Since powder samples have usually been 
found to show relatively high levels of paramagnetic im- 
purities and/or defects, this comparison suggests that 
the larger spin gaps in the powder samples may be a 
real effect arising from termination of the spin chains 
by defects. One would indeed expect finite segments of 
alternating-exchange chains to exhibit larger spin gaps 
than for the infinite chain. 

We have also listed in Table | the intradimer ex- 
change constant determined from y(T) (Ref. ^), 31 P 
NMR Knight shift, 31 K(T) (Ref. ||), and inelastic 
neutron scatteringo measurements of polycrystalline 
samples of the d 1 S =1/2 vanadium dimer compound 
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FIG. 2. Detailed crystal structure of AP-(VO) 2 P20 7 as 
viewed along the b-axis (top panel) and along the a axis (bot- 
tom panel). The rectangular boxes are outlines of the unit cell 
in the respective planes. The large spheres are V atoms and 
the small spheres are P atoms. Oxygen atoms (not shown) 
are at the vertices of the VO5 square pyramids and the PO4 
tetrahedra. The pairs of V atoms coupled by exchange con- 
stants Ji and J2 along the S = 1/2 AF alternating-exchange 
Heisenberg chains are shown. 



(VO)HPC>4 • 5H2O (vanadyl hydrogen phosphate hemi- 
hydrate), which is a precursor for the synthesis*, of, 
and has structural similarities to, AP-(VO)2P207.Q'EI In 
particular, the neutron scattering study of this com- 
pound confirmed the importance ami strength of the 
V-O-P-O-V superexchange pathway.E2l 

From the above neutron scattering measurements on 
single crystals of AP-(VO)2P20y, a second spin gap at 
a larger energy of 67 K was found in addition to the gap 
of 36 K_for coherent one-magnon propagation along the 
c axis. til The two spin gaps cannot both arise from one- 
magnon excitations in an isolated alternating-exchange 
chain and and the larger one was suggested to arise from 
neutron scattering from two-magnon triplet bound states 
of such chains, although the scattered neutron inten- 
sity was, larger than expected from scattering from such 
states.tLil A 2D model incorporating both nonfrustrating 
(J a , see Fig. |l|) and diagonal frustrating (J x ) AF inter- 
actions between alternating-exchange chains was subse- 
quently proposed.Lj The intrachain and interchain ex- 



^ c HP-(VO) 2 P 2 7 





FIG. 3. Detailed crystal structure of HP-(VO) 2 P 2 7 . The 
designations are the same as in Fig. ^. 



change constants were determined by a (very good) fit 
to the one-magnon dispersion relation (including the 
one-magnon spin gap of 36 K) measured in the inelas- 
tic neutron scattering measurements and by using the 
resultal3 of high-T x(^) measurements. The two alters, 
nating AF intrachain exchange constants so determinedEfl 
are listed in Table ffl. The two AF interchain interactions 
were predicted to, be rather large: J a /Ji = 0.203 and 
Jx/Ji — 0.255.C3 However, the x(T") calculated using 
this model in Ref. was found to be in poor agreement 
with experimental single-crystal x(^) dataE3 even at the 
highest measured temperature T ~ 200 K, which is sig- 
nificantly above the temperature (~ 70 K) at which x{T) 
shows a broad maximum and where the prediction would 
be expected to be quite accurate. On the other hand, we 
find in Sec. 1IIB 2 below that their x(T) prediction is in 
quite good agreement with experimental data if the com- 
parison is done in a somewhat different way. Additional 
calculations indicated that the frustrating AF interchain 
interaction stabilizes tworiaagnon bound states and thus 
supported the conjecturetil that the higher energy mode 
at 67 K is a triplet two-magnon bound state.O The same 
2D model containing alternating-exchange chains and 
frustrating AF interchain couplings and also a model con- 
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taining ferromagnetic interchain couplings were studied 
in Ref. 25 where the former model was found to agree 
better with the experimental neutron scattering data. 
From a Raman scattering study of single crystals, the 
spin-phonon interaction was suggested to be responsible 
for the formation of the two-magnon triplet bound states 
identified in the neutron scattering experiments, ratber 
than arising from frustrating interchain interactions.Ej 

Electron spin resonance (ESR) intensity versus tem- 
perature I(T) data obtained for a single crystal were in- 
terpreted in terms of an isolated AF spin dimer model, 
yielding a spin gap of 67 K,E3 corresponding to the higher- 
energy gap seen in the inelastic neutron scattering data. 
In this ESR study, transitions between the Zeeman lev- 
els within the one-magnon triplet band were reportedly 
not observed. These I(T) data were subsequently re- 
interpreted within the above 2D coupled alternating- 
exchange chain model as arising from transitions between 
the Zeeman levels within the one-magnon band and good 
agreement between the theoretical prediction for I(T) 
and the experimental I(T) data was found.Efl 

On the other hand, recent 31 P and 51 V NMR and 
magnetization versus applied magnetic field M(H) mea- 
surements at high fields and low temperatures have in- 
dicated that there are two magnetically distinct types 
of alternating-exchange V chains in AP-(VO)2P207, in- 
terpenetrating with each other, each with its own spin 
gapEScJ The two spin gaps inferred for the two-types 
of chains, 35(2) K from 31 K(T) measurementscJ and 
68(2) K from 51 V K(T) measurements^ and-33(l) K and 
62(3) K from M(H) measurements at 1.3K,c3 agree well 
with the above two spin gaps found from the neutron 
scattering measurements, respectively, thus providing an 
alternate explanation for the larger of the two spin gaps. 
Using additional information from the neutron .scattering 
one-magnon dispersion relation measurements, til the two 
alternating, exchange constants in each chain have been 
estimatedc3 and are listed in Tabic |. This model is not 
supported by the Raman scattering results.E9 However, 
recent unpublished inelastic neutron scattering measure- 
ments on a large single crystal at T = 2 K in zero and high 
magnetic fields showJJiat the two types of magnetic exci- 
tations found earliertil split in a magnetic field according 
to expectation for two triplet bands, results which con- 
stitute independent evidence for.-the validity of the two- 
chain model for AP-(VO)2P2C>7.c3 In addition, two pre- 
viously undetected magnetic excitations termed "shadow 
bands" were found which are thought to arise from the 
staggered alignment of successive V dimers along the V 
alternating-exchangc-chains within the structure (see the 
top panel of Fig. |).E§E3 

At present there is thus no universal agreement about a 
Hamiltonian which can self-consistently explain the var- 
ious experimental measurements probing the magnetism 
in AP-(VO)2P2C>7, although all of the models consid- 
ered recently contain 5 = 1/2 AF alternating-exchange 
Heisenberg chains as an essential element and a consen- 



sus is emerging that the two-chain model can explain 
many of the observed properties. 



B. HP-(VO) 2 P 2 7 

The high-pressure phase HP-(VO)2P20y was recently 
synthesized by heating polycrystalline AP-p(yO)2P207 
for 1/2 h at 700 °C under a pressure of 2 GPaE3 As shown 
in Fig. H and by comparison with Fig. ^, HP-(VO)2P20y 
has a simpler structure than AP-(VO)2P207. In 
HP-(VO)2P207, all V atoms are crystallographically 
equivalent but the same basic-structure as in the ambient- 
pressure phase was found.Ejl23 The similarities between 
the two structures suggest that HP-(VO)2P207 also con- 
tains 5 — 1/2 AF alternating-exchange chains, but of a 
single type.Ell Indeed, high-field M(H) measurements at 
1.3 K revealed ajSingle spin gap of ~ 23 K, consistent with 
this hypothesis Modeling of x(T) data below 30 K was 
carried out using the low-T approximation in Eq. (|7a|) be- 
low for the spin susceptibility of a ID spin system with a 
spin gap, yielding a similar spin gap of 27K.E3 The x(T) 
data above 30 K were analyzed using the spin susceptibil- 
ity of the 5=1/2 alternating-exchange Heisenberg chain 
model, yielding the exchange constants listed in Table | 
and the same spin gap A/fce = 27KE9 These estimates 
of the gap value are similar to the one-magnon gap of 
w 36 K in AP-(VO)2P207 found from the inelastic neu- 
tron scattering and other measurements discussed above. 

Equation ( |7a| ) has also been previously used to fit x(T) 
data for other 5 = 1/2 ID compounds, but to our knowl- 
edge all such studies, with the exception of a study of 
the Cu +2 r-i> = 1/2 two-leg Heisenberg spin ladders in 
SrCu203,Eil have assumed that A and the spin gap A are 
independently adjustable parameters when fitting exper- 
imental x(T) data. We discuss in Sec. [n] below that A 
is uniquely related to A and is not an independently a' 
justable parameter for a given type of ID spin lattice 
On the other hand, if one assumes that the spin lat- 
tice in a material has a spin gap but the type of ID 
spin lattice is unknown, or if the fit is not carried out 
in the low-temperature limit, then A and A would have 
to be treated as independently adjustable parameters. 
In the present work we evaluate the accuracy of using 
Eq. ( [7a] ) to determine the spin gap of the ID 5 = 1/2 
alternating-exchange Heisenberg chain from xC^O data 
for HP-(VO)2P207 by comparing the spin gap obtained 
using this approximation with the spin gap obtained by 
modeling the .same data set using the accurate theoreti- 
cal predictions for the spin susceptibility of the 5=1/2 
AF alternating-exchange Heisenberg chain model. 

Most recently, single crystals of HP-(VO)2P207 have 
been grownB3 The anisotropic spectroscopic splitting 
factors (g factors) of the 5 = 1/2 V +4 ions in a sin- 
gle crystal were determined from ESR measurements and 
the anisotropic piagnetic susceptibilites of a single crystal 
were measured.123 Here we perform detailed modeling of 
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these x{T) data using the g factors determined from ESR 
and using the 5 = 1/2 AF alternating-exchange Heisen- 
berg chain model for the spin susceptibility. We deter- 
mine the exchange constants within the chains and from 
these we obtain an estimate of the spin gap. An inde- 
pendent estimate of the spin gap is obtained by modeling 
only the low-temperature data. 



C. Plan of the Paper 

The plan for the remainder of this paper is as follows. 
In Sec. [Il| a summary is given of the theory that we will 
need to carry out and discuss the modeling described 
above. General considerations for fitting experimental 
data by theoretical predic tions for the spin susceptibility 
are discussed in Sec. Ill A[ New and literature x(T ) data 
for AP-(VO)2P207 are presented and fitted in Sec. IIIB. 
High- precision fits to the x( T) of a high-purity powder 



sample are presented in Sec. IIIB1. In this section, we 



show how the fitted exchange constants and spin gap(s) 
of the model of alternating-exchange chain(s) vary de- 
pending on whether a single-chain or two-chain model 
is used to fit the data, and on whether the g value is 
fixed or allowed to vary during the fits. The influences 
of interchain couplings on the exchange constants and 
spin gaps inferred from modeling x(T) for t he pow der 
sample are quantitatively determined in Sec. HI B 2 us- 
ing a molecular field theory for the interchain couplings, 
where we also compare our derived interchain couplings 
with the corresponding theoretical predictions of Uhrig 
and NormandEj which were obtained using a one-chain 
model. The x(T) data for two single crystals of AP- 
(VOWP2O7 are anal yzed using the two-chain model in 
Sec. pll B 3| . In Sec. [HIC we test our predicted disper- 
sion relations for the two proposed chains by comparison 
with the results of inelastic neutron scattering measure- 
ments. The x{T) d ata fo r HP-(VO)2P2 07 are presented 
and modeled in Sec. |lIIPj . The most accurate and precise 
x{T) data for this phase were obtained for a sample of 
crushed crystals. T hese dat a are analyzed using the one- 
chain model in Sec. 
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The influences of interchain 
couplings on the derived exchange constants and spin gap 
of this sample are considered in Sec. [Ill D 2 . Our evalua- 
tion of the accuracy of the spin gap obtained using a the- 
oretical low-T approximation to the spin susceptibility of 
a ID 5 = 1/2 spin system with a spin gap, as previously; 
used to analyze p owder x(T) data for HP-(VO)2P2C>7 
is given in Sec. HID 3. The anisotropic x(^) date 



for a single crystal of HP-(VO)2P207 are modeled in 
Sec. [II D 4. The powder average of the se single crys- 



tal x(^) data are modeled in Sec. HID 5 by two low-T 
approximations for the spin susceptibility to obtain an 
independent estimate of the spin gap. A summary of our 
modeling results and our conclusions are given in Sec. IV. 



II. THEORY 

The Hamiltonian for the 5=1/2 alternating-exchange 
Heisenberg chain is written in three equivalent ways as 

H = /J J\S2i-\ ■ S-2i + J 2 S 2 i • S2i+i (la) 

i 

= X] ^i^ai-i • S 2 i + aJiS 2l ■ S 2 i+i (lb) 

i 

= 5^ J(l + 5)S 2i -i ■ S 2t + J(l - 5)S 2t ■ S 2l+1 , (lc) 

i 

where 



Ji = J(l + S) = 



2J 

1 + a ' 



■h i-s 

a = — = 



J i 1 + 5 ' 

_ Ji _ _ Ji — J 2 _ I - a 
J ~ 2 J ~~ 1 + a ' 



J = 



■h 



(2a) 
(2b) 
(2c) 
(2d) 



with AF couplings Ji > J 2 > 0, < {a, S) < 1. The uni- 
form undimerized chain corresponds to a = 1, 5 = 
and J i — J 2 = J, whereas the isolated dimer has 
a = 0, 6=1, J 2 = and intradimer exchange inter- 
action J\. The form of the Hamiltonian in Eq. ([lc]), in 
which the appropriate variables are 5 and the average 
exchange constant along the chain J instead of a and 
the maximum exchange constant J\ as in Eq. ([Tb|), is of- 
ten used for compounds in which the spin dimcrization 
is weak and/or for systems showing a second-order spin 
dimerization transition with decreasing temperature such 
as occurs in spin-Peierls systems. 

The spin gap A for magnetic spin excitations from the 
5 = ground state to the lowest-lying 5 = 1 triplet ex- 
cited states for the alternating-exchange chain is uniquely 
related to the alternation parameter a and the larger ex- 
change constant J\ (or equivalently to S and J). The 
ratio A(a)/Ji for the 5=1/2 AF alternating-exchange 
Heisenberg chain was computed to high (< 1%) accu- 
racy for < a < 0.9, in a increments of 0.1, using-pulti- 
precision methods by Barnes, Riera and Tennant.Eia They 
found that their calculations could be parametrized very 
well by the simple expression 



A(a) 
•h 

or equivalently 



(1 - a) 3 / 4 (l + a) 1 ' 4 



J 



2 5 Z / A 



(3a) 



(3b) 



An expression for A(<5)/J which is thought to be more 
accurate (±0.0002) over the entire range < 5 < 1, ob- 
tained by fitting numerical A(S)/ J data by a generalized 
form of Eq. (pq) , isQ 



G 



J 



= 2S v{5) , 



(4a) 



where the exponent y(S) is given by 
y(S) = t/(l) + ni tanh 



In 5 ( In 5 
In 

rri! V m 2 



n 2 tanh" 





/ In S \ 






mi 


\m 2 ) 



(4b) 



with parameters 
2/(1) = 0.74922 , ni = 0.00776 , n 2 = -0.00685 



m 1 = 3.3297 , m 2 = -2.2114 



(4c) 



The expression for A(S)/ J in Eqs. (|J) can be transformed 
into an expression for A(a)/J± using the conversion ex- 
pressions (^). 

For notational convenience, we define the reduced spin 
susceptibility x* , reduced temperature t and reduced spin 
gap A* as 



X spil Vi 



t : 



Ji 



A* = — 



A 



(5) 



where x spin is the spin susceptibility, N is the num- 
ber of spins, fee is Boltzmann's constant and \* de- 
pends on both t and the alternating exchange parameter 
a = J 2 1 J i ■ 

An accurate but unwieldy two-dimensional function 
X*(t,a) for the S = 1/2 AF alternating-exchange 
Heisenberg chain has been derived for the entire range 
< a < 1 of the alternation parameter by a global 
fit to numerical quantum Monte Carlo (QMC) simula- 
tions and transfer-matrix density-matrix renormalizatiop. 
group (TMRG) and Bethe ansatz x*(t, a) calculations^ 
which we will not reproduce here but will explicitly use 
to model x(T) data for both AP- and HP-(VO) 2 P 2 7 . 
The absolute accuracy of this function for < a < 1 and 
0.01 < t is estimated to be < 2 x 10~ 4 , which corresponds 
to < 0.1% of the susceptibility at the broad maximum. 
For practical purposes of fitting experimental x(T) data, 
this function can be considered toJae exact for t > 0.01. 

Troyer, Tsunetsugu, and WiirtzO have derived a gen- 
eral expression for the low-T limit of x*(t) f° r a one- 
dimensional spin system with a spin gap, assuming that 
(i) the one-magnon dispersion relation is nondegencr- 
ate (apart from the Zeeman degeneracy), (ii) the lowest 
magnetic excited states are one-magnon S — 1 triplet 
excitations, (iii) k^T <C A and k^T <C one-magnon 
bandwidth, and (iv) the one-magnon dispersion relation 
E(k) is parabolic near the minimum according to (in the 
present notation) 



(6) 



where k is the wavevector in the direction of the ID sys- 
tem and a is the (average) nearest-neighbor spin-spin dis- 
tance. These assumptions hold for the present case of the 
5=1/2 AF alternating-exchange Heisenberg chain ex- 
cept for the limit_p — as discussed below and for a = 1 
for which-A = OB With these four assumptions, X*(t) is 
given bytfJ 



with 



-A*/t 



.4 



(t < A*, bandwidth/ Ji) (7a) 



(7b) 



The dimensionless dispersion parameter c* in Eq. (|g) 
has a unique relationship to the reduced spin gap A* for 
any given ID spin system. For example, according to the 
model of Ref. |^, this relationship for the S = 1/2 AF 
alternating-exchange Heisenberg chain gives the value of 
the parameter A in Eq. (J7a|) as 



A = 



(8) 



/2¥/(A*) ' 

which in turn yields the low-T limit of x* {t) in Eq. (|7a|) 



as 

X *(t) = ^=l (— ) 12 e- At/t , (t<A*) (9a) 

V ' V2^/(A*) V t J v ' v ' 

where the dimensionless function /(A*) is the solution of 

/2( A *) 



E 



2A* 



(9b) 



and E(x) is the complete elliptic integral of the second 
kind. 

From Eq. (||), the parameter A in Eq. ( |7a| ) is not 
an independently adjustable parameter but instead is 



was 



a unique function of the reduced spin gap A* 
also previously inferred for two-leg spin ladders. Eil In ad- 
dition, we see from Eq. (|9a|) that the two independent 
parameters of x*(t,ot) can be written for low tempera- 
tures as x*(V A *: A *) = X*(k B T/A, A/Ji). Finally, and 
perhaps most importantly, the high-temperature limit of 
the low-temperature regime in which x* V.is closely ap- 
proximated by Eqs. is of order A*/10.H At such low 
temperatures xCO is immeasurably small, and hence the 
spin gap obtained by analyzing experimental x{T) data 
for various compounds up to temperatures corresponding 
to a sizable fraction of A* using Eq. ((7a|) leads to fitted 
spin gap values which may be significantly different from 
the actual spin gaps. 

The result for the low-i limit of X*(t) in Eq. ( ffa| ) is 
not valid for the isolated dimer, which is one limit of 
the alternating-exchange chain with a = J 2 — 0, because 
the assumption (iii) that k B T <C one-magnon bandwidth 
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required for that equation to hold is violated at all finite 
temperatures. For the isolated dimer, the one-magnon 
bandwidth is identically zero and the reduced spin gap is 
A* = 1. The reduced spin susceptibility is given exactly 
by 



,*, dimer 



(<) 



i(3 + eV*) 



with the low-temperature limit 



X 



*,dii 



c (t) = \e- 1/t (t«l) 



(10a) 



(10b) 



The temperature depend ence of the prefactor to the ex- 
ponential term in Eq. ( |l0b|) is different from that in 
Eq. ( [fa] ) . As discussed in Ref . [|, a crossover occurs with 
decreasing temperature at low temperatures in the effec- 
tive prefactor from a 1/t dependence to a 1/yt depen- 
dence if < a <C 1. 

As noted above, the form for the low-T behavior of the 
spin susceptibility in Eq. (|7a|) is valid only at very low 
temperatures. Many years ago, Bulaevskii found that 
his numerical values of x*(*> a ) f° r the S = 1/2 AF 
alternating-exchange Heisenberg chain, computed from 
an analytic theory based on the Hartree-Fock approxima- 
tion, could be fitted over a-*elatively large temperature 
interval 0.033 < t < 1/4 byB 



x*(t) 



t 



(11) 



and he tablulated A and A* versus the alternation pa- 
rameter a. A recent extensive numerical studyO of his 
theory confirmed that the numerical predictions of his 
theory in the above-cited low-temperature range are fit- 
ted better by the form (|i~i"| ) than by ((7a|). In addi- 
tion, this study showed that although the fitting param- 
eter A* (a) approximately follows the actual spin gap of 
Bulaevskii's theory, significant discrepancies occurB Fi- 
nally, a detailed numerical comparison of the prediction 
of Bulaevskii's theory for x*(t> a ) with QMC simulations 
and TMRG calculations of this quantity showed that 
Bulaevskii's theory is unsuitable for accurately extract] 
ing a values from experimental x(T) data when a < l.u 
For Heisenberg spin lattices consisting of identical spin 
subsystems with susceptibility Xo(t) which are weakly 
coupled to each other, the molecular field theory (MFT) 
prediction for the reduced spin susceptibility x* (f) m t he 
paramagnetic state of the system is 



**(*) = 



Xo(t) 



l + AxS(<) 



or equivalently 



1 



+ A, 



X*(t) Xo(*) 
where the MFT coupling constant A is given by 



(12a) 



(12b) 



A = 



/ -j Jmax 
3 



(12c) 



the prime on the sum over j signifies that the sum is only 
taken over exchange bonds from a given spin Si to 
spins Sj not in the same spin subsystem, and J max is the 
exchange constant in the system with the largest magni- 
tude. By definition, the expression for Xo(*) does not 
contain any of these Jij interactions which are external 
to a subsystem. Within MFT, Eqs. ( |i~2| ) are correct at 
each temperature in the paramagnetic state not only for 
bipartite AF spin systems, but also for any system con- 
taining subsystems coupled together by any set of FM 
and/or AF Heisenberg exchange interactions. 



III. MODELING OF EXPERIMENTAL X (T) DATA 

A. Introduction 

We fitted the x(T) data per mole of V spins-1/2 in 
(VO)2P207 by the general expression 



X(T) = X o + 



a 



imp 



T-6 



+ X spin (T) 



imp 



with 



and 



Xo = X 



vv 



(13a) 



(13b) 



X spin (T) 



(0.3751 



cm 3 K 



(k B T\ 

mol I J\ j fee ^ ^ J\ ' 



V 



X 



(13c) 



where N\ is Avogadro's number, /ie is the Bohr mag- 
neton, fee is Boltzmann's constant and g is the spectro- 
scopic splitting factor (g factor) appropriate to a partic- 
ular direction of the applied magnetic field with respect 
to the crystal axes. 



The first term xo in Eq. (13a), according to Eq. ( |13b| ), 
is the sum of the nearly isotropic orbital diamagnetic 
atomic core contribution ^ corc and the anisotropic or- 
bital paramagnetic Van Vleck contribution x VV which 
are normally nearly independent of T. Using the values 



-12 and -47 ^0" 6 cm 3 /mol for V+ 4 , O 



and (PO4) 3 , respectively,! 3 ^ we obtain 



= -6.1 x 10 



-5 



cm 
mol V 



(14) 



for (VO)2P2C>7. The second term in Eq. (13a) is an ex- 
trinsic impurity Curie- Weiss term with impurity Curie 
constant C lmp and Weiss temperature 9i mp which gives 
rise to a low-temperature upturn in x(^) which is not 
predicted by theory for the third term, the intrinsic spin 
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susceptibility x spin (T), and is assumed to arise from finite 
chain segments containing an odd number of spins, im- 
purity phase intergrowths in the crystals, paramagnetic 
impurity phases and/or defects. The Ci mp and #i mp pa- 
rameters can be anisotropic if the paramagnetic impurity 
principal directions are fixed with respect to the crystal 
axes, as can occur in a singl e crystal of a material such as 
studied here in Sec. HID 4, rather than being randomly 
distributed. 

Unless otherwise stated, we assume that the spin sus- 
ceptibility x spnl ( T) in Eq. (13a), written in terms of 
X*(t,a) in Eq. ( |l3c| ), is the intrinsic spin susceptibil- 
ity per mole of spins- 1/2 in an AF alternating-exchange 
Heisenberg chain. The explicit expr essi on for x*(t, a) of 
this chain is given in Ref. [5| In Eq. (13c), Ji is the larger 
of the two ( J i and J 2 with J\ > J 2 > 0) AF alternating 
exchange constants along the alternating-exchange chain, 
as denoted in Sec. || above. 

One of the parameters entering the calculated spin 
susceptibility 



'CO 



Eq. 



(13c) is the g value 
of the V magnetic moments. Measurements of the 
anisotropic g values of the spins in both AP-(VO)2P2C>7 
and HP-(VO)|gP|a£)7 have been carried out using ESR 
measurementalJ'ES and the results are listed in Table ||. 
The significant differences between the g values of the 
two phases of (VO)2P2C>7 reflect the differences in the 
local bonding of the V atoms with the coordinating O 
atoms in the two structures. Comparison of the average 
g value for HP-(VO)2P207 with those for V in the "txcllis. 
layer" compounds i?V 2 5 with R = Ca, Mg and Na,ETt3 
also shown in Table Q, suggests that the local crystalline 
electric field (CEF) at the V sites in HP-(VO) 2 P 2 7 is 
closer to that in these compounds than to the CEF in 
AP-(VO) 2 P 2 7 . 

A measure of the goodness of a fit to experimental 
x{T) data is the statistical x 2 P er degree of freedom 



> 
o 



AP-(VO) 2 P 2 7 Powder 
H = 1 kG 




50 100 150 200 250 300 350 
T (K) 

FIG. 4. Magnetic susceptibility \ versus temperature T 
for a powder sample of AP-(VO) 2 P 2 7 (o) from 2 K to 350 K. 
The solid c urve is a two-dimensional fit to the 357 data points 
using Eq. (L3a) assuming a spin susceptibility x spm (T) con ~ 
sisting of those of equal numbers of two independent isolated 
5 = 1/2 antiferromagnetic alternating-exchange Heisenberg 
chains. The dashed curve is the fitted x spi "(T). The fit- 
ted exchange constants and the derived spin gaps for the two 
chains are listed in the column labeled "Fit 6" in Table H . 



DOF 



1 Np 

7pEW T f)-»H(Ti)] ! 



N n -P 



(15) 



where N p is the number of data points in the data set 
and P is the number of independent fitting parameters. 
This is the quantity that is minimized during our non- 
linear least-squares fits to experimental x(T) data. An 
additional measure of the quality of a fit is the relative 
rms deviation <7 rms of the fit from the data, given by 



TABLE II. g-factors parallel (<?||) and perpendicular 
(g±) to the principal local crystalline electric field and/or 
crystal structure axis and the powder-averaged value 
g — \J (gf + g\ + C/D/3 for V +4 S = 1/2 species in several 
vanadium oxide compounds. Samples are polycrystalline un- 
less otherwise noted. The literature references are given in 
the last column. 



Compound 


511 


9± 


g 


Ref. 


AP-(VO) 2 P 2 7 


1.94 


1.98 


1.97 




17 


crystal 


1.937(a) 


1.984(6, c) 


1.969 




22 


HP-(VO) 2 P 2 7 


1.928(l)(a) 


1.974(1)(6), 


1.958 




32 


(crystal) 




1.971(l)(c) 








CaV 2 5 






1.957(1) 




37 


MgV 2 5 






1.96 




38 


NaV 2 5 (crystal) 


1.938(2) 


1.972(2) 


1.961(2) 




3£ 


crystal 


1.936(2) 










crystal 


1.95 


1.97 


1.96 




y 


crystal 


1.936 


1.974, 1.977 


1.962 




42 











1 p 

— y 



x( T i) - Xfit(Ti)" 1 



x(Ti] 



(16) 



The fits were all carried out on a 400 MHz Macintosh G3 
computer using the software Mathematica3.0. 



B. Magnetic Susceptibility of AP-(VO) 2 P 2 7 



1. Powder sample 

The purpose of the present section is to test con- 
sistency with experimental x(T) data of the aaQdel of 
Yamauchi and co-workers for AP-(VO) 2 P207,Ea'c3 dis- 
cussed in the Introduction, in which this compound is 
proposed to consist magnetically of equal numbers of two 
independent types of isolated S — 1/2 AF alternating- 
exchange Heisenberg chains with spin gaps of about 35 
and 68 K, respectively. 
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The x{T) of a polycrystalline ("powder") sample of 
AP-(VO)2P207 of mass 172.2 mg and with a moss-green 
color was measured from 2 to 350 K in a magnetic field 
of 1 kG and the results are shown as the open circles 
in Fig. ^. The details of the sample preparation will 
be presented elsewhere. The color of the sample indi- 
cates that it is stoichiometric with a vanadium oxida- 
tion state very close to +4.E3 The sample was a cylinder 
of 4 mm diameter and 7 mm length. There was no dif- 
ference between field-cooled and zero- field-cooled x(T) 
measurements. The quality of the sample, judging from 
the very small Curie- Weiss upturn at low temperatures 
due to magnetic impurities and/or defects, is better than 
previously reported for any powder sample and is about 
the same as recently reported for a high-quality single 
crysta£3 as shown in Figs. |^ and ||(b) below. 

We will describe in detail our modeling results for this 
sample to indicate how the parameters and the quality of 
fit change for various types of fits. Similar variations were 
found from modeling x(T) data for two AP-(VO)2P207 
single crystals below. To make contact with previous 
modeling of x(T) for this material, we first fitted the data 
by Eqs. ( p"3| ) assuming that x spm (^) is due to a single type 
of S = 1/2 AF alternating-exchange Heisenberg chain, 
where the g value is either fixed at the powder-averaged 
value g — 1.969 determined from ESR measurements for 
AP-(VO) 2 P 2 7 as shown in Table g yielding "Fit 1", 
or allowed to vary during the fit ( "Fit 2" ) . Throughout 



the modeling in this section, we use the expression (Q) 
to determine the spin gap from the fitted exchange con- 
stants for an alternating-exchange chain. The parameters 
obtained from each fit are shown in Tabic III, together 
with the statistical % 2 /DOF and a lms for each fit. The 
defect and/or impurity Curie constant is equivalent to 
the contribution of about 1.6mol % with respect to V of 
spins 1/2 with g — 2. The values of the alternating 
exchange constants J\ and J2 and the spin gap A are 
respectively similar in each fit and are about the same 
as previously estimated from similar fits to x(T) data for 
this compound, respectively (see Table |). Next, we fitted 
the x(^) data using the above r«j|©del of Yamauchi and 
co-workers for AP-(VO)2P207,ESE3 but where we con- 
strained the spin gaps to be 36.2 and 66.7 K as found, 
from the inelastic neutron scattering measurements, til 
and again either fixed g = 1.969 ("Fit 3") or allowed 
g to vary during the fit ("Fit 4"). In order to enforce 
the constraint on the two spin gaps using the expression 
for A(i5)/J in Eqs. (Q), it was more convenient to use J 
and (5 as the independent parameters in x spm (^) during 
our least-squares fits instead of J\ and a. The param- 



eters obtained from the two fits are shown in Table III 



together with other parameters derived from the fitted 
ones. As can be seen from the values of the x 2 /DOF and 
<7 rms , the qualities of the two fits show dramatic improve- 
ments over those of Fits 1 and 2 where only a single type 
of alternating-exchange chain was assumed in the model- 



TABLE III. Fitted and derived parameters for x{T) of a high-quality powder sample of AP-(VO)2P207. A derived quantity 
is marked by an asterisk (*). A quantity with a "ee" in front of it was constrained to have the value listed and was not fitted. 
Quantities and are quantities associated with two independent isolated chains, respectively. The g value in Fits 1, 

3, and 5 were constrained to be the powder-averaged value 1.969 from ESR measurements (see Table [n]), whereas in Fits 2, 

4, and 6 the g value was fitted. Fits 1 and 2 assume a single type of isolated alternating-exchange chain, whereas Fits 3-6 
assume the presence of two independent isolated alternating-exchange chains. In Fits 3 and 4 the two respective spin gaps were 
constrained to have the values found (Ref. [ill ) from inelastic neutron scattering measurements, whereas in Fits 5 and 6 the two 
spin gaps were each allowed to vary independently. Our favored fit parameters are those of Fit 6. 



Quantity 


Fit 1 


Fit 2 


Fit 3 


Fit 4 


Fit 5 


Fit 6 


Xo (10- 5 ^) 


-1.8(1) 


-4.5(6) 


-1.37(6) 


-5.2(1) 


-1.99(2) 


-2.3(1) 


CW (io- 3 ^f) 


5.8(1) 


6.3(2) 


3.72(2) 


4.52(3) 


4.26(2) 


4.29(2) 


imp (K) 


-6.4(2) 


-6.9(3) 


-3.08(4) 


-3.95(3) 


-3.89(2) 


-3.91(2) 


9 


= 1.969 


1.993(6) 


= 1.969 


2.006(1) 


ee 1.969 


1.972(1) 


4 1] /kB (K) 


130.9(1) 


132.4(4) 


137.0(8)* 


144.4(4)* 


128.7(4) 


130.3(7) 


jfVfe (K) 






122.7(6)* 


126.0(4)* 


129.9(4) 


128.8(5) 


J«/fcB (K) 


97.9(2)* 


99.8(5)* 


118.2(6)* 


121.7(3)* 


107.6(5)* 


109.5(9)* 


4 2) /k B (K) 






76.4(5)* 


75.2(2)* 


83.7(4)* 


82.6(6)* 


Q (D 


0.748(1) 


0.754(2) 


0.863(1)* 


0.8728(5)* 


0.836(1) 


0.840(2) 


Q (2) 






0.622(2)* 


0.6178(7)* 


0.644(1) 


0.641(2) 


J (1) /fe (K) 


114.4(2)* 


116.1(5)* 


127.6(6)* 


136.2(4)* 


118.1(4)* 


119.9(8)* 


J (2) Ab (K) 






99.6(5)* 


98.5(2)* 


106.8(4)* 


105.7(6)* 


5 m 


0.1442(6)* 


0.140(1)* 


0.0735(5) 


0.0679(3) 


0.0893(6)* 


0.087(1)* 


5 (2) 






0.233(1) 


0.2363(5) 


0.2165(7)* 


0.219(2)* 


A«/fe (K) 


53.6(3)* 


53.2(4)* 


ee 36.2 


ee 36.2 


38.8(3)* 


38.6(5)* 


A^/k B (K) 






ee 66.7 


ee 66.7 


67.7(4)* 


67.5(5)* 


X 2 f 10 -5 cm 3 ^ 2 
DOF V moivy 


1.31 


1.24 


0.194 


0.0379 


0.0130 


0.0117 


CTrms (%) 


1.97 


1.94 


0.617 


0.226 


0.138 


0.132 
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ing. However, the g value obtained from Fit 4 is some- 
what larger than expected. 

Finally, we fitted the same data set in Fig. [| us- 
ing the above mod el of Yamauchi and co-workers for 
AP-(V0)2P207,c3O where we again either fixed g — 
1.969 ("Fit 5") or allowed g to vary during the fit 
("Fit 6"), but where we did not constrain the fitting 
parameters J\ and a of the two independent chains to 
yield the respective spin gaps found from the inelastic 
neutron scattering measurements. The fitted parameters 
are listed for each fit in Table III, together with the sta- 
tistical x 2 /DOF and cr rms and the derived A for each 
fit. We checked that the identical fitted parameters are 
obtained for Fit 5 independent of whether the starting 
parameters are the fitted parameters of Fit 1, for which 
the exchange constants in the two chains are identical, 
or of Fit 3 for which they are different. Fit 6 is shown as 
the solid curve in Fig. || and the fitted x spm (T) is shown 
as the dashed curve. 

The fitted and derived parameters for Fits 5 and 6 in 
Table [II exhibit a number of important features. First, 



the qualities of Fits 5 and 6 to the data are far superior 
to those of Fits 1 and 2. Second, the values of the al- 
ternation parameters and spin gaps for the two indepen- 
dent isolated chains of the model did not converge to the 
same respective values for the two chains, but rather are 
clearly differentiated. Third, the fitted g value from Fit 6 
is identical within the respective errors with the powder 
averaged g value in Table || determined from ESR mea- 
surements. Fourth, and perhaps most importantly, the 
two spin gaps derived from the respective exchange con- 
stants for the two chains are respectively nearly identical 
to the two spin gaps found from the single-crystal inclas= 
tic neutron scattering measurements by Garrett et al.t3 
and with the values inferred previously by Yamauchi and 
co-workers from high-field magnetization measurements 
and a subset of the NMR measurements. EjcH It seems 
very unlikely that the two spin gaps we deduce from this 
model could be so close to those determined from other 
independent measurements without the model being es- 
sentially correct. The exchange constants and spin gaps 
we derived from x(^) data for the same powder sample 
in H = 50 kG, data which are not otherwise discussed 
here, are identical within the respective errors to those 
we obtained above for H = 1 kG. Finally, we will see 



in Sec HI D 1 below that when the two-chain model is 



used to extract the exchange constants within the pro- 
posed single-chain high-pressure phase HP-(VO)2P207, 
essentially the same exchange constants and spin gaps are 
obtained for both chains of the model. This result indi- 
cates that our fitting procedure can clearly differentiate 
between pairs of chains that have the same or different 
spin gaps, respectively. 

We conclude that our analysis of x{T) is precisely 
consistent with the model of Yamauchi and co-workers 
for the nature of the important spin interactions in 
AP-(VO)2P2C>7. Our values of the spin gaps of the 
two independent isolated alternating-exchange chains of 



the model are in good agreement with those deter- 
mined from thfuc. high-field magnetization and NMR 
measurementscj'EH and with the two values determined 
from the inelastic neutron scattering measurements, til re- 
spectively. 



2. MFT analysis of interchain coupling 

The neutron scattering measurements on single crys- 
tals showed unambiguously that interchain coupling J a 
along the a axis of the structure, perpendicular to the 
alternating-exchange chains and parallel to the legs of the. 
structural two-leg ladders (see Fig. |l|), is not negligible. til 
However the ratio \J a /Ji\ was estimated from fits to the 
data to be only 2-3%, where J\ is the larger of the two 
exchange couplings alongjthe alternating-exchange chains 
running along the c axis. til Another estimate can be ob- 
tained as one-half the ratio of the average total disper- 
sion of the two presumed one-magnon bands in the a axis 
direction [16(4) K] to the one-magnon excitation energy 
along the direction of the alternating-exchange chains 
at the zone boundary (180 K), yielding a slightly larger 
\J a /Ji\ ~ 4.5(10)%. This interchain coupling along the a 
axis was of course ignored in the fits to the experimental 
x(T) in the previous section. Here we obtain an estimate 
of the strength of this interchain coupling J a from anal- 
ysis of the powder x(^) data presented in the previous 
section. In the absence of accurate calculations of x* (t) 
for this case, we will utilize the prediction of MFT given 
in Eqs. ( |l2|) for the influence of the interchain coupling 
on x{T). 

In order to apply Eqs. (|2|) to the present modeling 
framework in which two distinct types of alternating- 
exchange chains (1) and (2) are assumed to be present, 
one must appropriately define the "isolated subsystem" 
discussed in Sec. ||. Here, an isolated subsystem con- 
sists of one set of the two chains (1) and (2). Thus, the 
reduced susceptibility of our isolated subsystem is 



Xo 



Ng^i 
Jf 8 



(i) 



where 



IJ1 



J (1) ' 



-^chain 



f(2) 



(;«««) + 



Till 



-n 



(2) 



^cha 



2) nX2) 



■n 



(2) 



Jl 



max 



t(1) t(2) 

,j 1 , ,j 1 



(17) 



Then the MFT cou pling constant A in Eqs. ( |l2| ) is, ac- 
cording to Eq. (12c), the average interchain coupling in 
the a axis direction of a spin in one of the two dis- 
tinct alternating-exchange chains with all spins in the 
respective adjacent alternating-exchange chains. Assum- 
ing that a spin in each chain is coupled to two nearest 
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FIG. 5. Parameters Jj 1 ' 2 ' and a' 1,2 ^ and rms deviation 
cr rms of the fits to the powder x{ T ) data for AP-(VO) 2 P20 7 
in Fig. Q by MFT for coupled alternating-exchange chains of 
types (1) and (2) versus the MFT interchain coupling con- 
stant A. A g = 1.969 was assumed in all of the fits. The spin 
gaps A*- 1 ' 2 ' for the two distinct chains (1) and (2), derived 
from the respective J^ 1 ' 2 ' and a' 1,2 ' values, are also plotted. 
The lines connecting the data points are guides to the eye. 



rameters J^ 1 ' 2 ' and a^ 1 ' 2 ) for the chains (1) and (2), 
the rms deviation cr rms of the fit from the data, and 
the spin gaps A' 1 ' 2 ) for the chains (1) and (2) de- 
rived using Eqs. (Q), are plotted versus A in Fig. || for 
—0.4 < A < 0.4 in A increments of 0.05, where positive 
(negative) values of A correspond to AF (ferromagnetic 
FM) coupling J a . Not plotted in Fig. || are the fitted 
values of xo, Cj mp , and 9- lm p, which for A = —0.4, 
and 0.4 are -7.0(1), -1.99(2), 3.8(2) x 10" 5 cm 3 /mol V, 
4.65(7), 4.26(1), 3.4(1) xl0~ 3 cm 3 K/molV, and -4.0(1), 
-3.89(2), -3.0(2) K, respectively. From Fig. §, a pro- 
nounced minimum (0.095 %) occurs in a rms at A « —0.05, 
which corresponds according to Eq. ( |l8| ) to a FM J a 

with \J a /J^'\ ps 0.025. This is-quantitatively consistent 
with the above-cited estimatestiJ of this ratio based on 
a one-chain model with anisotropic (in spin space) spin 
interactions for the one-magnon dispersion relations ob- 
served by inelastic neutron scattering in Ref. |ll|. The 
fit to the data for A = —0.05 is shown as the solid 
curve m Fig. H where xo = -2.50(2) x 10~ 5 cm 3 /mol V, 
Cimp = 0.00428(1) cm 3 K/molV, 6> imp = -3.87(2) K, 
jW/fc B = 139.5(2) K, a« 



0.8597(4), j[ 2) /k B = 
123.7(2) K, = 0.6319(6), A^/k B = 37.5(2) K, and 
A( 2 '//cb = 66.1(2)K. In deriving the spin gap for each 
chain using Eqs. (0) we have implicitly assumed that the 
spin gap is unaffected by the interchain couplings. 

A somewhat more precise estimate of A is obtained by 
allowing this parameter to vary during the fit. The fit 
parameters and derived spin gaps of the two chains (1) 
and (2) for the best fit are 



Xo 



-2.37(4) x 10" 



mol V 



Cimp = 0.00427(1) 



cm 3 K 
mol V ' 



np = -3.87(2) K , A = -0.037(4) , 



7 (i) 

—3— = 137.1(7) K , = 0.855(2) 
Kb 



(19a) 



r(2) 

-J— = 124.8(4) K , a (2) = 0.634(1) , 
Kb 



DOF 



0.89 10 



molV 



0.099% , 



neighbor spins in the a direction by exchange constant 
J a , one obtains 



A = 2- 



Ja 



JTD 



(18) 



We fitted the data in Fig. ^ by Eqs. (13), assum- 
ing g — 1.969 for all of the fits and using Eqs. (12a), 
(pd), and fQ, where t = min(£ (1) , t (2) ), to determine 
the spin susceptibility x spin (T). The resulting fitted pa- 



A« 
Kb 



37.8(5) K , 



A( 2 ) 
Kb 



66.4(3) K . 



The spin gaps are similar to and may be compared with 
those for Fit 5 in Table III for the same data, in which g — 
1.969 was also assumed but where A = 0. From Eq. ( p^ ) 
which assumes a nearest neighbor interchain coordination 
number of 2, we obtain the average interchain coupling 
strength 
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X = 0.916, One-Chain Model 



a = 5.6 % 

rms 




AP-(VO) 2 P 2 O y 
Powder, H = 1 kG 







50 100 150 200 250 300 
T (K) 

FIG. 6. Magnetic susceptibility \ versus temperature T for 
a powder sample of AP-(VO)2P207 (°) from Fig ^. The solid 
curve is a fit to the data using Eq. (13a) assuming a spin 
susceptibility x spm (r) consisting of those of equal numbers of 
two 5 = 1/2 antiferromagnetic alternating-exchange Heisen- 
berg chains which are coupled using MFT with a ferromag- 
netic coupling constant A = —0.05. The dotted curve is the 
MFT prediction using the exchange constants found for the 
one-chain model by Uhrig and Normand (Ref. pl| ), for which 
the MFT coupling constant is strongly antiferromagnetic with 
the value A = 0.916. 



Ja 

fc R 



(1) 



-2.5K 



(19b) 



As noted in the Introduction, Uhrig and NormandEfl 
proposed a model for AP-(VO)2P207 in which only one 
type of alternating-exchange chain occurs and in which 
the (AF) interchain couplings are given by 



Ji 



124K, a = 0.793, 



Ja 
Jl 



0.203, 



Jx 
J: 



0.255 , 
(20a) 



where these parameters were obtained by fitting the 
one-magnon inelastic neutron scattering dispersion rela- 
tion datallil for tha-Iower band and using the observed 
Weiss temperaturetj obtained by fitting experimental 
x(T) data at high temperatures by a Curie- Weiss law 
(as predicted by MFT). Assuming that the interchain 
couplings do not affect the spin gap, Eqs. (0) yield 



A 

— = 44.0 K . 
Kb 



(20b) 



This i s ab out 16% larger than our estimate for A' 1 ) in 
Eqs. (19a). Of course, since their exchange constants 
were determined by fitting their theory to the experi- 
mental neutron scattering data, their spin ga p is t he ob- 
served value (w 36 K) and not that in Eq. (20b). The 
discrepancy arises because they find that the spin gap 
does depend o n the interchain couplings, as further dis- 
cussed in Sec. Ill C below. 



Since the interchain spin coordination number for each 
of the interchain couplings is 2, the value of the MFT 
interchain coupling constant predicted by Eq. (12c) is 



Ji Ji 



(20c) 



Inserting the parameters of Uhrig and Normand in 
Eq. ( |20a| ) into this equation yields A = 0.916. Our fit 
parameters and their variations with A in Fig. ^ argue 
against this very large AF value of A. To further illus- 
trate the discrepancy within MFT between this one-chain 
theory and the experimental x(T) data, shown as the dot- 
ted curve in Fig. is the predicted x(T) using g = 1.969, 



the 



Ao 



and 



values obtained for A 



-0.05 in 



Sec. IIIB2, and the MFT prediction for the spin suscep- 
tibility in Eqs. jl^), where Xo(*) is that of the isolated 
alternating-exchange chain for which the exchange con- 
stants estimated by Uhrig and Normand in Eq. (20a) 
were used. The relative deviation of the prediction from 
the data is cr rms = 5.6 %, which is about 60 times larger 
than the cr rms obtained using the two-chain model for 
A = —0.05. The agreement of both theoretical predic- 
tions with the data at high temperatures is expected and 
in fact is required for either model, since the MFT is 
most accurate at high temperatures where it yields the 
Curie- Weiss law. The significant differences between the 
predictions of the two models only become apparent at 
the lower temperatures. 

In summary, our high-precision fits to the xCO data 
using the model of two independent chains, in which 
nearest neighbor chains along the a axis are coupled us- 
ing MFT, indicate that the average interchain coupling 
is weakly ferromagnetic, in agreement with the analysis 
of neutron scattering data by Garrett et alc3 using a 
one-chain model and in disagreement with the one-chain 
model of Uhrig and NormancO with strong AF inter- 
chain couplings. 



3. Single crystals 

In this section we analyze the anisotropic x{T) data 
for two single crystals of AP-(VO)2P207. The data 
for a small dark green crystal of mass 2.0(1) mg ("crys- 
tal 1" ), measured in a magnetic field of 5 T using a Quan- 
tum Design SQUID magnetometer, have not been re- 
ported previously. There was no discernable difference 
between field-cooled and zero-field-cooled x(T) data for 
this crystal. The data for crystal 2 were reported by 
Prokofiev et al. in Fig. 4(a) of Ref. p2|and were measured 
in a magnetic field of 2 T. An overview of the anisotropic 
x(T) data for the two crystals is shown in Fig. [?], where 
the x(T) of the powder sample in Fig. [| is shown for 
comparison as the dashed curve. The data for the two 
crystals are in agreement on a coarse scale. The powder 
averages of the data for both crystals lie above the data 
for the powder sample for T > 25 K, although in the case 
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FIG. 7. Overview of the anisotropic magnetic susceptibil- 
ities x versus temperature T for crystals 1 (open symbols) 
and 2 (filled symbols, from Ref. of AP-(VO) 2 P 2 7 along 
the a axis (circles), b axis (squares) and c axis (diamonds). 
Also shown for comparison is x(T) for the powder sample 
from Fig. kl (dotted curve). 
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FIG. 8. Fits to the anisotropic magnetic susceptibilities x 
versus temperature T for crystals 1 (a) and 2 (Ref. ^2|) (b) of 
AP-(VO)2P207. The fits are shown as solid curves and the 
spin susceptibilities as dashed curves. Note that the temper- 
ature scales in (a) and (b) are different. 



of crystal 1 this difference is not significant since the 
uncertainty in the crystal mass is about 5 %. Crys- 
tal 1 shows a larger Curie- Weiss-type paramagnetic de- 
fect and/or impurity contribution at low temperatures 
than crystal 2. Quantitative differences are seen between 
the anisotropic x(T) for the two crystals. In particu- 
lar, above about 30 K the Xa (T) of crystal 1 agrees with 
Xb,c(T) of crystal 2 and Xb,c{T) of crystal 1 agrees with 
Xa(T) of crystal 2. These qualitative anisotropy differ- 
ences cannot arise from inaccuracy in, e.g., the crystal 
masses, which would only affect the respective ordinate 
scale. 

The a-, b- and c-axis x(T) data for each crystal in Fig. |?] 
were fitted simultaneously using Eqs. ( |i~3| ) by writing x a s 
a diagonal tensor. We assumed the two-chain model for 
the spin susceptibility of each crystal where the values 
of J\ and a for each chain are the same for all three 
crystal directions, the anisotropic g values are the same 
for both chains, and allowed xo, Cimp and 6*j mp to be 
different for each chain and for the three field directions, 
for a total of 16 fitting parameters. The 4D fits obtained 
for crystals 1 and 2 are shown as the sets of three solid 
curves in Figs. |](a) and ||(b), respectively, and the fitted 

TABLE IV. Fitted and derived parameters for the 
anisotropic x(T) °f our crystal 1 and for crystal 2 (Ref. |%]) of 
AP-(VO)2P207. A derived quantity is marked by an aster- 
isk (*). Quantities and are the quantities associated 
with two independent isolated chains, respectively. 



Quantity 


Crystal 1 


Crystal 2 


XOa (10- 5 ^) 


-3.4(3) 


-5(3) 


Xoo (10- 5 ^) 


-1.6(2) 


0(2) 


Xoc (10- 5 ^) 


-7.7(2) 


1(2) 


r fin" 3 ™i 3 K\ 
lmpa mo iy) 


14.9(1) 


7.1(5) 


Cimpf, (10 molv ) 


14.7(1) 


3.7(4) 


n. (in -3 cm3K ) 

lmp c \w molv ) 


13.8(1) 


7.7(5) 


Oimpa (K) 


-9.85(8) 


-10.9(9) 


#impb (K) 


-10.06(8) 


-4.5(5) 


#impc (K) 


-9.73(8) 


-10.7(6) 


9a 


2.003(2) 


1.98(2) 


9h 


1.984(2) 


2.00(2) 


9c 


1.946(2) 


1.98(2) 


4 1} /k B (K) 


122(3) 


125(14) 


j[ 2) /k B (K) 


143(4) 


128(11) 


a m 


0.803(9) 


0.80(5) 




0.648(6) 


0.65(4) 


4 1} /k B (K) 


98(4)* 


101(16)* 


4 2) /kB (K) 


93(3)* 


83(12)* 


j (1) Ab (k) 


110(3)* 


113(15)* 


j (2) Ab (k) 


118(3)* 


105(12)* 


s m 


0.109(5)* 


0.11(3)* 


s m 


0.214(5)* 


0.21(2)* 


A (1) //cb (K) 


42(3)* 


43(11)* 


A^/k B (K) 


74(3)* 


66(10)* 


X 2 ( 1Q -5 cm 3 \ 2 
DOF V molV/ 


0.0384 


2.0 


(Trms (%) 


0.159 


1.61 
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parameters for both crystals are given in Table [V where 
the goodnesses of fit for the two crystals are also listed. 
The anisotropic spin susceptibilities )d£ m (T) were de- 



rived using Eq. ( 13a ), i.e., by subtracting the respective 
Xo and defect and/or impurity Curie- Weiss terms from 
the %(T) fit function, and are plotted versus temperature 
for crystals 1 and 2 as the two sets of three dashed curves 
in Figs. ||(a) and ||(b), respectively. From Eqs. (|l3|), the 
only source of anisotropy in ^ s P m (T) is the anisotropy in 
the g factor. Also listed in Table [TV] are the spin gaps 
computed using Eqs. (^) for the two distinct alternating- 
exchange spin chains in each crystal. 

Several features of the data in Table are of note. 
First, as qualitatively expected from Fig^TJ the concen- 
tration of paramagnetic defects and/or impurities in 
crystal 1 is about a factor of two larger than in crystal 2. 
Second, the spin gaps of the two chains in each of the 
crystals 1 and 2 are consistent within the error bars with 
each other and with those found in the above section 
for the high-quali ty p owder sample of AP-(VO)2P207 



as listed in Table III 



The large error bars on the fit- 
ted parameters for crystal 2 arise-in large part because 
the resolution in x f°r the dataE3 above 20 K is only 
1 x 10 -5 cm 3 /mol V, which corresponds to a relative res- 
olution of, e.g., 1% at 20 K and 0.5% at 70 K. Third, 
the fitted g a values are similar to, but differ in detail 
from, the corresponding ESR values for AP-(VO)2P207 
in Table ||. These discrepancies between the respective 
g a values may originate at least in part from the large 
Curie- Weiss contribution in crystal 1 and the low resolu- 
tion of the x data for crystal 2. 
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FIG. 9. Comparison of our predicted dispersion relations 
of the proposed two types of alternating-exchange chains 
in AP-(VO)2P207 with the dispersion relations along the 
chain direction measured using inelastic neutron scattering at 
T = 10 K by Garrett et al. (•, Ref. [H]). The two solid curves 
are the dispersion relations predicted from our exchange con- 
stants determined from fits to x{T) of a high-purity powder 
sample by the isolated chain model, and the dashed curves are 
the corresponding curves for a MFT-coupled chain model. In 
each case, the dispersion relations were calculated from our 
intrachain exchange constants using the Fourier series to tenth 
order in a given by Knetter and Uhrig (Ref. ^). The dotted 
curve for each chain is the dispersion relation in Eqs. ( |2^ ) 
incorporating the interchain coupling according to the model 
of Uhrig and Normand (Ref. |^). 



C. Dispersion Relations of AP-(VO) 2 P 2 7 

The one-magnon dispersion relation E{k c ) / J\ in the 
chain direction for the isolated S — 1/2 AF alternating- 
exchange Heisenberg chain, with and without a frustrat- 
ing second-neighbor coupling, was recently-calculated to 
tenth order in a by Knetter and Uhrig. c3 Thus it is 
possible to make a direct and accurate comparison of 
the dispersion relations predicted for the two proposed 
alternating-exchange chains on the basis of our exchange 
constants in AP-(VO)2P207 with those determined by 
Garrett et a/.tLil using inelastic neutron scattering (INS) 
measurements. For this comparison, we first use the ex- 
change constants for the two chains from Fit 6 in Table III 
determined from our fit to the x{T) data for the high- 
purity powder sample using the isolated chain model. 
The predicted dispersion relations for the two chains are 
shown as the two solid curves in Fig. ^ where the ex- 
perimental INS data (•) are also plotted. Also shown as 
dashed curves are the dispersion relations predicted for 
the two cha ins using the intrachain exchange constants 
in Eq. (19a) found from our fit to the same x(T) data 
by the same two-chain model but where the chains are 
coupled along the a axis using MFT. The range of our 



prediction for the dispersion relation of each chain is thus 
approximately given by the region between the respective 
pair of solid and dashed curves. 

Our predicted dispersion relations for the two chains 
agree well with the experimental inelastic neutron scat- 
tering data in Fig. ^ except near the zone boundary at 
k c = tt/c. Our prediction is that two peaks in the scat- 
tered neutron intensity versus energy at this wavevec- 
tor should be seen with an energy separation of about 
20-40 K (w 2-3 meV) , contrary to the single data point 
at this wavevector in Fig. ^|. However, the error bar on 
the data point in Fig. at k c = tt/c is only the accu- 
racy of determining the position of the centroid of the 
neutron scattering peak-and is not a direct measure of 
the width of the peak.Ej After most of the fits to the 
x(T) data described in this paper and the determinations 
of the exchange constants were completed, we learned 
that recent unpublished INS measurements on a large 
crystal of AP-(VO)2P207 indeed show two peaks at this 
wavevector with an energy splitting of about 2 meV,c3 
which partially confirms our prediction. We also predict 
from Fig. |^ that the dispersion relations of the two chains 
should cross within the intermediate wavevector regime. 
We are not aware of experimental INS data that address 
this aspect of the dispersion relations. 
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FIG. 10. Fits to the dispersion relations of the proposed 
two types of alternating-exchange chains in AP-(VO)2P207 
in the a direction and perpendicular to the chain direction, 
measured using inelastic neutron scattering at T = 10 K by 
Garrett et al. (•, Ref. [H]). The two solid curves are fits to 
the respective data by Eqs. (^l|), which is the special case for 
k c — of the general dispersion relation calculated by Uhrig 
and Normand (Ref. |^). 

We note that the two sets of exchange constants in 
Table IV for the two chains in crystal 1 also predict an 
energy splitting of the neutron peak at k c — n/c, but 
in this case the predicted dispersion relations of the two 
chains do not cross. We consider the experimental data 
and modeling for the high-purity powder sample to be 
more r eliable than for this crystal, for reasons discussed 
in Sec. Ill B 3 above, and hence expect that the dispersion 
relations of the two chains will ultimately be observed to 
cross. 



A quantitative estimate of the interchain couplings J a 
and J x can be obtained by fitting the observed^ disper- 
sion relation for each of the two chains in the a direc- 
tion, shown in Fig. [l0| by the prediction of Uhrig and 
NormandE-3 

E(k a ,k c = 0) = A(Ji,a,/i+,//_) 

oo 

+ J\ a n (a, /Lt+, fi-)[cos(nk a a) — 1] (21a) 

n=l 

where fi± = (J a ± J x )/Ji. To third order in a, 
and/or ji + their general dispersion relation yields 



«i 



/'- 



4 + a(2 + ^+)-a 2 - fx 2 _ 



8 



a 3 



(21b) 



We found that the data in Fig. |l0] could not be fitted 
by Eqs. (|2l|) assuming J x = and J a /Ji ~ —0.019 as 



inferred from A ~ —0.037 in Eqs. (|l9|). However, we can 
still retain this experimentally determined value of A by 
allowing J x to be nonzero. In particular, from Eq. (20c) 
one obtains /i + = A/2. Therefore we set /i+ = —0.019 
and used the experimental ly d etermined values of J\ and 
a for each chain in Eq. (19a) in the fit to the respec- 
tive dispersion relation. This left A and /i_ as the only 
adjustable parameters. We did not use the experimen- 
tally determined A values in Eqs. (|l9|) because the total 
width of the experimental dispersion relation for each 
chain is relatively small and the small difference between 
the experimental A value and the neutron scattering re- 
sult would give a significant systematic shift to the fit 
(see also below). The fits are shown as the solid curves 
in Fig. |To| for chains (1) and (2), respectively, for which 
the parameters are 



(i) 



Ji 



t(2) 



-0.035 , 



= -0.032 



J 



(i) 



Ji 



-0.016 



A« 



J 



(2) 



= +0.013 



A( 2 ) 



36. 2K , 

= 65. 7K . 

(22) 



Thus we find that a small but finite AF value of the inter- 
chain interaction J x is necessary to fit the dispersion data 
perpendicular to the alternating-exchange chains if we 
retain the fitted A value in Eqs. (|l9|). Since the nearest- 
neighbor interchain interaction J a is ferromagnetic, the 
next-nearest-neighbor AF interchain interaction J x is not 
a geometrically frustrating interaction. 

On the other hand, an equally good and nearly identi- 
cal fit for each chain as shown in Fig. [l(] can be obtained 
assuming that J x = if we relax the above condition 
on A. In this case we still use the exchange constants 
in Eqs. (|l9| ) but we set /i_ = /i + = fi in the fit to the 
transverse dispersion data for each chain, yielding the 
same respective gap values as in Eqs. (|22|) , but where 
= -0.050(2) and ^ = -0.044(2), so that the 
Ja/ Ji = M an( l A = 2/i w —0.10 values are larger in 
magnitude than in Eqs. (|22]). 

The dispersion relation par 
lated by Uhrig and Normand, 

interchain couplings are included to third order in a, /i 
and/or yields 



el to a chain ( 1 1 c) calcu- 
in which the J a and J x 



E(k a = 0, k c ) = A + [E (k a = 0, fe c ) - B o (0, 0)] 

oo 

+ Ji b n (a, /Lt+, fi-)[cos(nk c c) — 1] (23s 



with 

h -- 



fi-a 



a 
1 



3/i_a 2 
16 



(23b) 



where Eq(Q, k c ) is the dispersion relation for the, isolated 
alternating-exchange chain to tenth order in a.LJ The re- 
sulting dispersion relations calculated for the two chains 
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using Eqs. (|23|), with the the intrachain exchange con- 
stants and spin gaps in Eqs. ( |l9| ) and the interchain cou- 
plings given in Eqs. (22), are shown as the dotted curves 
in Fig. ^. These are respectively very similar to those 
for the MFT-coupled chain parameters already plotted 
as the dashed curves in Fig. g. 

An inconsistency in our fit to the transverse dispersion 
relation for each chain is that the calculated spin gap 
is smaller than the ob serv ed and fitted one in Fig. |lC| . 
The spin gap in Eq. (21a) is given by the third order 
expansion of Uhrig and NormandE.il as 



HP-(VO) 2 P 2 ? Crushed Crystals 



A = Ji- 



A (a) [i- 



4 + a(2 + /i+) - a 2 



Ji 4 

2 3 

+ ^-(4-2a + ^ + )-^ 



(24) 



where A (a)/Ji is the spin gap in Eqs. (^) or (Q) in the 
absence of interchain coupling. For the exchange con- 
stants in Eqs. ( |19a| ) and (22) used to fit the transverse 
dispersion relations in Fig. 1C for chains (1) and (2), from 
Eq. (||) we obtain AW/fe = 29.2 K and A^/k B = 
59.8 K, respectively. These spin gaps are each signifi- 
cantly smaller than the fitted ones in Eqs. ((H), respec- 
tively. These discrepancies arise because the interchain 
couplings change the spin gap, contrary to our implicit 
assumption when we fitted the experimental x(^) data 
using the MFT-coupled chain model, so the J\ and a 
intrachain exchange parameters for each chain derived 
from the MFT fit to these data must be considered in 
the present model to be effective values. The degree to 
which the effective exchange constants differ from the ac- 
tual ones is difficult to evaluate. The combined analysis 
we have done of the susceptibility and dispersion rela- 
tions is as rigorous as can be done without having in 
hand an accurate theoretical expression for the spin sus- 
ceptibility of the system which includes the influence of 
interchain couplings and concommitant changes in the 
two spin gaps. 



D. Magnetic Susceptibility of HP-(VO) 2 P 2 7 

1. Crushed crystals 

We begin our modeling of x(T) for HP-(VO) 2 P 2 7 
using data shown in Fig. [if] which we obtained for a 
72.2 mg sample of crushed green transparent single crys- 
tals. These data are expected to be more accurate and 
yield more reliable values of the exchange constants and 
spin gap than the fits to the data for a powder and 
for a very small single crystal discussed in the follow- 
ing two sections, respectively. The data were modeled 
using Eqs. ( |i~3| ) in which x*(t) is the theoretical reduced 
susceptibility for the S — 1/2 AF alternatine=exchangc 
Heisenberg chain as proposed by Azuma et alEQ The fit 
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FIG. 11. Magnetic susceptibility \ versus temperature T 
for a sample of HP-(VO)2P2C>7 crushed single crystals (o) 
from 2 to 400 K. The solid curve is a two-dimensional fit to 



the 300 data points using Eq. (13a) and the S — 1/2 antifer- 
romagnetic alternating-exchange Heisenberg chain model for 
the intrinsic spin susceptibility x spa> (T)- The dotted curve is 
the fitted x spm (r). The fitted values of the g factor, the larger 
of the two exchange constants in the alternating-exchange 
chain Ji, the alternation parameter a = J2/J1, and the de- 
rived spin gap A are listed. 



TABLE V. Fitted and derived parameters for %{T") of 
crushed single crystals of HP-(VO)2P207 obtained using 
the one-chain and two-chain models. A derived quantity is 
marked by an asterisk (*). Quantities A™ and are the 
quantities associated with two independent isolated chains, 



respectively. 


Quantity 


One- Chain Model 


Two-Chain Model 


Xo (10- 5 ^) 


-0.1(4) 


-1.0(3) 


« m p (10- 3 


4.88(6) 


4.70(4) 


0imp (K) 


-1.99(5) 


-1.79(3) 


9 


1.951(3) 


1.958(3) 


j[ r> /k B (K) 


135.2(3) 


134(14) 


JpVfe (K) 




136(14) 


a m 


0.8727(4) 


0.90(2) 






0.85(2) 


J«/feB (K) 


118.0(3)* 


120(15)* 


4 2) Ab (k) 




115(15)* 


j (1) Ab (k) 


126.6(3)* 


127(15)* 


j (2) Ab (k) 




126(14)* 


s m 


0.0680(3)* 


0.054(10)* 


5 (2) 




0.082(10)* 


A (1) /fc B (K) 


33.9(2)* 


29(7)* 


A^/kB (K) 




39(7)* 


X 2 f 10 -5 cm 3 ^ 2 
DOF V moivy 


0.36 


0.16 


<T me (%) 


0.58 


0.35 
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is shown as the solid curve in Fig. 11 and the fitted 
X spm (T) is shown as the dotted curve. The fitted g, 
Ji, and a values are listed in the figure, along with the 



MFT Fits to Crushed Crystal x(T) data 



The other parameters 
The fitted g value is 



spin gap computed using Eqs. ( 
of the fit are shown in Table 

very close to the powder average value 1.958 in Table |l| 
obtained from ESR measurements. The impurity Curie 
constant is equivalent to the contribution of 1.3mol% 
with respect to V of spins 1/2 with g — 2. 

We used the data set in Fig. [ll| to estimate typical non- 
statistical errors that may arise when using the two-chain 
model to fit the x(T) data for AP-(VO) 2 P 2 7 in the 
above two sections. In addition, if the present one-chain 
model for HP-(VO)2P20y is appropriate, then a fit by 
the two-chain model should yield very similar exchange 
constants and spin gaps for the two chains of the model, 
which ideally would be respectively identical for the two 
chains. The parameters of the two-chain fit are com- 
pared with those of the above single-chain fit in Table [v|. 
We see that the fitted parameters of the two chains us- 
ing the two-chain model are the same within the limits 
of error with each other and with the parameters of the 
single-chain model, respectively. This result indirectly 
confirms that the large differences between the exchange 
constants and spin gaps found above for the two chains 
in AP-(VO) 2 P 2 7 are reliable. 



2. MFT analysis of interchain coupling 

There are no data available for the strength of the 
interchain coupling in HP-(VO)2P207. W e estim ate 
this coupling in the same way as in Sec. [II B 2| for 
AP-(VO)2P207. The most precise and accurate x(T) 
data available for HP-(VO)2P207 are those for the 
crushed crystal sample presented and discussed in the 
previous section. We determined the fitting parameters 



Xo, Cimp, #imp, Ji, and a as a function of the MFT 
interchain coupling constant A over the range —0.5 < 
A < 0.5, where the fixed g = 1.958 was assumed in 
all these one-chain model fits. The cr rms and the Ji 
and a parameters are plotted versus A in Fig. 12, along 



with the spin gap A determined from J\ and a using 
Eqs. (||). At A = -0.5, 0, and 0.5, the fitted param- 
eters xo, Cimp, and lmp were, respectively, —6.1(2), 
-0.86(9), 3.7(2)xlO- 5 cm 3 /molV, 5.52(9), 4.96(4), 
4.53(8) xlO" 3 cm 3 K/molV, and -2.33(9), -2.03(4), 
— 1.81(8) K. From the top panel of Fig. [l^, the cr rms 
shows an approximately parabolic variation with A, with 
a minimum at A k —0.05, indicating a weak ferromag- 
netic interchain coupling as also deduced above for AP- 
(VO) 2 P 2 7 . 

We next allowed A to vary during the fit to determine a 
more precise value. The fit parameters and derived spin 
gap of the alternating-exchange chain for the best fit are 




-0.4 



-0.2 



0.0 

X 



0.2 



0.4 



FIG. 12. Parameters Ji and a and the rms deviation 



<j rms of the fits to the y(T) data for crushed crystals of 
Fig. 



HP-(VO) 2 P2 07 in Fig. |llj by MFT for coupled alternat- 
ing-exchange chains versus the MFT interchain coupling con- 
stant A. A g = 1.958 was assumed in all of the fits. The spin 
gap A, derived from the Ji and a values for each A, is also 
plotted versus A. The error bars are shown for Ji, a, and A. 
The lines connecting the data points are guides to the eye. 
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Xo = -1-4(2) x 10" 



molV 



Cimp = 0.00502(4) 



cm 3 K 
mol V ' 



np = -2.06(4) K , A = -0.055(14) , 



^- = 136.2(3) K , a = 0.8748(5) , 

fc B 



(25a) 



DOF 



a T1 



0.59% , 



kB 



= 33.8(2) K 



The spin gap is identical with that obtained for the one- 
chain fit in Table [v] for the same data, in which g = 1.958 
was also assumed but where A = 0, and the other param- 
eters are also very similar, respectively. From Eq. (18), 
we obtain the average interchain coupling strength 



Ja 



AJi 

2fcB~ 



-3.7 K . 



(25b) 



3. Powder sample: low-T fits 

The x(T) of a powder sample of JJP-(VO)2P207 was 
previously reported by Azuma et al.ffj shown as the open 
circles in Fig. [l3]. A fit of the data up to 30 K by 
Eqs. (|l3]), where X*(t) is the low-T approximation for the 
spin susceptibility of a gapped ID S = 1/2 spha, system 
in Eq. ( (7a[ ) , yielded the spin gap A/k B = 27 KB In this 
fit, the parameter A in Eq. ( [7a]) was treated as an inde- 
pendently adjustable parameter. In this section we carry 
out a precise fit of the same data set by Eqs. (|l3|) using 
the accurately known x*{t,a) spin susceptibility predic- 
tion far the 5 = 1/2 AF alternating-exchange Heisenberg 
chain The present fit was carried out in order to com- 
pare the fitted parameters respectively obtained from the 
two types of fits to the same x(T) data set for the same 
sample. 

We fitted the x(T) data in Fig. [jl using Eqs. Q, 
where x*(^j Q! ) 1S that for the S = 1/2 AF alternating- 
exchange Heisenberg chain given in Ref. |5|, and with the 
g value fixed at the spherically- averaged value 1.958 de- 
termined from single-crystal anisotropic ESR measure- 
ments (see Table D). The resulting fit is shown as the 
solid curve in Fig7]l3|, and the fitted spin susceptibility 
is shown as the dashed curve. The parameters of the fit 
are 



Xo = 3.8(3) x 10 



-5 



cm 



Cimp = 0.0121(2) 



cm 3 K 
mol V 



molV ' 



= -2.6(1) K, 
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FIG. 13. Magnetic susceptibility \ versus temperature T 
for a powder sample of HP-(VO) 2 P20 7 (o) from 2 K to 400 K 
(Ref. ^) . The solid curve is a two-dimensional fit to the 164 
data points using Eq. (13a) and the S = 1/2 antiferromag- 
netic alternating-exchange Heisenberg chain model for the in- 
trinsic spin susceptibility X spm (T), using the powder-averaged 
g-value of 1.958 determined from ESR measurements. The 
dashed curve is the fitted X spm (T). The fitted values of 
the larger of the two exchange constants in the alternat- 
ing-exchange chain Ji, the alternation parameter a = J2/J1 
and the derived spin gap A are listed. 



-i = 141.6(4) K , a = 0.849(2) 



(26) 



The Cimp value in Eq. ( J26| ) is rather large, equivalent to 
the contribution of 3.2mol% of spins-1/2 with respect to 
V and with g = 2. Therefore, the fitted prefactor 1/ J\ to 
X*(t) m Eq. (13c) could be too small by about 3% if the 
magnetic species in the impurity phase is S = 1/2 V +4 . 
On the other hand, if the impurities/defects have a spin 
larger than 1/2, the influence on the fitted parameters 
could be much smaller. In order to test this possible 
influence, we next allowed the g value to be an adjustable 
parameter in the fit, which has the effect of allowing the 
amount of V in the HP-(VO)2P207 phase relative to that 
in the impurity phase to be variable. The parameters 
obtained, 



Xo = 0.0(12) x 10" 



Gm P = 0.0125(3) 



cm 3 K 
mol V 



cm 
mol V 



= -2.7(1) K, 



3 = 2.00(1), -i = 143.9(8) K 

Kn 



0.854(2) , (27) 



are very similar to those obtained in the above fit with 
fixed g. The fitt ed v alue of the prefactor 5 ,2 /(Ji/fce) 
to X*(t) in Eq. (13c) increased, as anticipated, from 
2.71 K" 1 to 2.77 K^T Although the fitted g value in- 
creased slightly from the value used in the first fit, it is 
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FIG. 14. (a) Magnetic susceptibility x versus temperature 
T for single crystal HP-(VO)2P207 along the a-axis (o), 6-axis 
(•) and c-axis (open squares) directions from 2 K to 300 K. 
(b) Expanded plot of the data in (a) at low temperatures 
below the broad maximum in x(T) at about 75 K. In (a) 
and (b), the set of three solid curves is a four-dimensional 
fit to all 885 d ata points for the a-, b- and c-axis direc- 
tions using Eq. (13a) and the S = 1/2 antiferromagnetic al- 



ternating-exchange Heisenberg chain model for the intrinsic 
spin susceptibility x spm (T). The dashed curves are the fitted 
Xa" n (T) for the magnetic field along the a = a-axis (short 
dash) and b- and c-axes (longer dash); the fitted Xb Pm (T) and 
Xc pm (T) are indistinguishable on the scale of the figure. 



still close to 2. We conclude that the magnetic impurities 
and/or defects giving rise to the Curie- Weiss term have 
little influence on the fitted exchange constants in the 
HP-(VO)2P207 phase in the sample. 

Taking J x /k B = 142 K and a = 0.85 from Eq. ((H), 
Eq. (pah yields the spin gap A/fee = 40 K, about 



50% larger than the above value of 27 K estimated by 
Azuma et al!Ei by fitting the same data up 30 K (i.e., up 
to a T ~ A/&b) using the low-T approximation to the 
spin susceptibility in Eq. (|7a|). Thus from the respective 
values of A, we find that the error arising from estimat- 



ing the gap value by fitting low-T x(T) data using the 
low-T approximation for the spin susceptibility is about 
50% in this case. The temperature range over which the 
low-T approximation is fitted to the exper imenta l data 
is expected to influence this error (see Sec. [IID5). 



4- Single crystal 

The x(T) data reported by Saito et aZ.H for a 0.26 mg 
single crystal of HP-(VO) 2 P20 7 are plotted in Fig. |l|. 
By comparison of these data with those for the powder 
sample in Fig. [l3|, the Curie- Weiss impurity and/or defect 
contribution to x(T) for the crystal at low temperatures 
is seen to be significantly smaller than for the powder 
sample. From our fit below, we find that it is in fact 
about a factor of three smaller. This much smaller impu- 
rity contribution enhances the reliability of the fitted ex- 
change constants and the derived spin gap obtained from 
modeling the data for the single crystal. At higher tem- 
peratures, the magnitude of the powder-averaged x(T) 
for the single crystal is very similar to the x(T) for the 
powder, as would have been expected. 

The modeling of the single-crystal x(T) data was car- 
ried out in a similar way as for the two crystals of 



AP-(VO)2P207 in Sec. HI B 3 , except that here we use a 
one-chain model instead of a two-chain model for the spin 
susceptibility. All 885 data points of the a-, 6-, and c-axis 
X(T) data sets in Fig. |l4| were fitted simultaneously using 
Eqs. (13) by writing x as a diagonal tensor and using the 
three fixed g values determined for fields along the Lhrcc 
principal axis directions from ESR measurements^ as 
given in Table |fj above, respectively. With eleven fitting 
parameters, the data to parameter ratio is 80. The four- 
dimensional fit obtained is shown as the set of three solid 
curves in Fig. [l4], and the fitted parameters are given in 

TABLE VI. Fitted parameters for the anisotropic mag- 
netic susceptibility of single crystal HP-(VO)2P207 along 
the a-, b-, and c-axis directions from 2 to 300 K using the 
one-chain model for the spin susceptibility. Derived quanti- 
ties are shown with an asterisk (*). The values of Ji, J2, a, 
J, S, and A are the same for all three crystal axis directions. 



Quantity 


a-axis 


b- axis 


c-axis 


xo(io- 5 ^rv) 


3.4(2) 


-0.9(2) 


0.0(2) 




3.40(6) 


3.90(6) 


4.39(7) 


0imp (K) 


-2.48(9) 


-1.98(8) 


-2.75(9) 


Ji/k B (K) 




131.6(1) 








0.8709(5) 




x vv (io- 5 ^) 


9.5* 


5.2* 


6.1* 


J 2 /k B (K) 




114.6(2)* 




J Iks (K) 




123.1(1)* 




8 




0.0690(3)* 




A/k B (K) 

X 2 ( W -5 cm 3 \ 2 
DOF V moivy 




33.4(2)* 






2.1 




(Time (%) 




1.43 
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Table [Vl] along with the goodnesses of fit. The spin sus- 
ceptibilities for the three crystal directions are plotted 
versus temperature in Fig. |j as dashed curves. 

The average of the three fitted Ci mp values is 3.9 x 
1CP 3 cm 3 K/mol V, which is equivalent to the contribu- 
tion of 1.0mol% with respect to V of paramagnetic 
species with S — 1/2 and g = 2. This contribution is 
about a factor of three smaller than that found for the 
powder sample in the previous section. The average of 
the three #i mp values is —2. OK, about the same as for 
the powder sample. The negative sign of 6>i mp may indi- 
cate AF interactions between the defect and/or impurity 
magnetic moments. The #i mp can also arise from single- 
impurity-ion CEF effects, and/or as a reflection of partial 
T-dependent paramagnetic saturation of the paramag- 



netic impurities at low temperatures in the fixed field of 
the measureme nts. The anisotropic x vv values are de- 
rived from Eqs. (|l3b | ) a nd ( |l4] ) using the fitted anisotropic 
Xo values in Tabl e |Vl| . The resulting three x vv values 
are listed in Table \V% 

From the above global fit to the anisotropic x(T) data 
for all three field directions, we find Ji/fes = 131.6(1) K 
and a = 0.8709(5). The average of the two exchange 
constants J± and Ji according to Eq. ( |2d| ) is J/fce = 
123.1(1) K, and the alternation parameter expressed in 
the form of 8 is obtained using our Ji and a values and 
Eq. (^c]) as 8 = 0.0690(3). Using our fitted a and J\ pa- 
rameters, Eq. ( |3a| ) yields the spin gap A/fee = 33.2(2) K. 
Using the derived J and 8 parameters, the more accurate 
Eqs. (§) predict the spin gap to be A/fc B = 33.4(2) K, 
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FIG. 15. Parameters obtained from fitting the powder-average x{T) data from Fig. [L4|for single crystal HP-(VO)2P207 from 
2 K up to a temperature T max by Eqs. ( |l3| ) using the expressions for the low-T spin susceptibility x*M i n Eqs. ( f7a| ) (o, "Fit 1") 
and (JfTI) (•, "Fit 2"), respectively. 
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which is the same to within the error bars as the first 
estimate. 

A summary of all the fitted and derived quantities ob- 
tained in this section from our modeling of the xCO mea- 
surements for single crystal HP-(VO)2P207 is given in 
Table 0. 



5. Powder average of single crystal x(T): low-T fits 

Additional fits to the powder-averaged single crystal 
x(T) data were carried out at low temperatures to ob- 
tain an estimate of the spin gap which is independent of 
the model for the gapped spin susceptibility. The powder 
average was used in order to reduce the number of pa- 
rameters needed to fit the data. We used the general fit 
expressions ( |l3|) in which the spin susceptibility X*W is 
given by the low-T approximations in Eq. ( f7a| ) ("Fit 1") 
or Eq. (|ll]) ("Fit 2"), and where the prefactor parameter 
A was fitted independently of the spin gap. The only 
difference between Fits 1 and 2 is the exponent of 1/t in 
the prefactor to the exponential in the expression for the 
low-T x*(t), which is 1/2 for Fit 1 and 1 for Fit 2. 

The powder-averaged single-crystal x(^) data from 
Fig. [l4| were fitted from 2 K up to a maximum tempera- 
ture T max , and the fitted parameters obtained from Fits 1 
and 2 are plotted versus T max in Fig. |l5| as open and filled 
circles, respectively. Also shown in Fig. [l5]are the statis- 
tical variances (x 2 per degree of freedom) obtained from 
both types of fits. The variances for both fits are quite 
similar and both have a minimum for T max w 25 K. How- 
ever, the impurity Curie constant and Weiss temperature 
for Fit 1 in Fig. 15 are strongly and nonmonotonically 
dependent on T n 
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HP-(VO) 2 P 2 7 Crystal: Powder Average 



in contrast to the corresponding de- 
pendences for Fit 2. In addition, the values of the fitted 
Xo values for Fit 1 in Fig. [l5| are all strongly negative. 
Since xo cannot be more negative than ^ corc as estimated 
above in Eq. (fl4|), because the \ VV in Eq. (13b) is neces- 
sarily positive, the Fit 1 fits for all the fitted T max values 
are unphysical and hence the other parameters obtained 
using Fit 1 are most likely also highly inaccurate. 

Shown in Fig. [l6] are the two optimum fits obtained 
for T max = 25 K for Fits 1 and 2, respectively, along 
with the respective fitted spin susceptibilities x Bpm (T)- 
The x spm (^') f° r the optimum Fit 1 is highly unlikely, as 
are the fit parameters as just noted. On the other hand, 
X spm (T) and the fit parameters for the optimum Fit 2 
are reasonable. The values of the fitted parameters for 
the optimum Fit 2 with T max = 25 K in Fig. |l| are 

Xo = 0(4) x 10" 



C imp = 0.0038(4) 



cm 3 K 



molV 



molV ' 



= -1.8(3) K, 



A = 0.118(2) 
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FIG. 16. Powder-averaged magnetic susceptibility x ver ~ 
sus temperature T below 50 K from Fig. [L4| for single crystal 
HP-(VO) 2 P20 7 (o). The solid and dotted curves are Fits 1 
and 2 to the data from 2 to 25 K using the low-T approxima- 
tions for the spin susceptibility in Eqs. (^) and ([j~l|), respec- 
tively. The corresponding fitted spin susceptibilities x spm (T) 
are shown as short- and long-dashed curves, respectively. Ex- 
trapolations of the curves to lower and higher temperatures 
are also shown. 

The first three parameters are very close to the corre- 
sponding powder-averaged anisotropic single crystal val- 
ues in Table VI which were obtained in the above section 
assuming the S = 1/2 AF alternating-chain Heisenberg 
model for the spin susceptibility, and the two spin gaps 
are nearly identical. The agreement between the spin 
gaps found from the two independent fits supports the 
applicability of this spin Hamiltonian to the spin system 
in HP-(VO) 2 P 2 7 . 



IV. SUMMARY AND CONCLUSIONS 



We have carried out detailed modeling studies of the 
magnetic susceptibilities x(^) of powder and single crys- 
tal samples of the ambient- and high-pressure phases of 
(VO)2P2C>7. The major goal of the modeling was to 
determine whether the recent proposals of a two-chain 
model for AP-(VO) 2 P 2 7 (Refs. fffi and |j) and a single- 
chain model for HP-(VO)2P207 (Ref. |30|) are consistent 
with the respective experimental x(^) data. Using the 
high-accuracy theoretical x*{t> a ) function for isolated 
5 = 1/2 AF alternating-exchange Heisenberg chains in 
Ref. §, high precision tests of these models were possible. 

The x{T) data for each phase were first analyzed us- 
ing the AF alternating-exchange chain model for isolated 
chains. We found that the proposed models are strongly 
supported by our high-precision fits to the x(T) data 
for each phase, from which the exchange constants and 
spin gap of each type of chain in each phase were de- 
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termined. We then considered the case of coupled chains. 
The influences of interchain couplings on the values of the 
intrachain exchange constants and the spin gap of each 
type of chain in the two phases were evaluated from addi- 
tional fits to the x(T) data where the interchain coupling 
was treated in the molecular field approximation. For 
both phases, we find that the interchain molecular field 
coupling constant is weakly ferromagnetic with a value 
A w —0.05, in agreement with Ref. [ll] and in disagree- 
ment with Ref. |24|. Assuming that the interchain coordi- 
nation number is two and using the values of A and the 
intrachain exchange constants, the interchain exchange 
coupling constant along the a axis direction is computed 
to be J a /k B w —3.0(5) K in both phases of (VO) 2 P 2 7 . 
Thus although our modeling of x(T), from which the in- 
trachain exchange constants and spin gaps were derived, 
did not explictly incorporate the influence of the magnon 
dispersion perpendicular to the chain direction, we be- 
lieve that our mean-field treatment effectively captures 
most of these effects on x(T) since the interchain cou- 
pling is found to be very weak compared to the intra- 
chain couplings. This supposition is confirmed by our fit 
to the low-T powder-averaged data for a single crystal of 
HP-(VO)2P207 by a model-independent low-T approx- 
imation for the spin susceptibility of a ID spin system, 
which yielded a spin gap of 33.8(2) K that is identical 
within the error bars to the spin gap of 33.4(2) K ob- 
tained from a fit to the three complete anisotropic x(T) 
data sets for the crystal using the x*(t) spin suscepti- 
bility function for the alternating-exchange chain. The 
good agreement of the respective spin gaps with those 
in Table [| obtained from inelastic neutron scattering and 
NMR measurements also supports the magnetic models 
that we have used for the two phases. 

According to the usual simple model for d orbitals of 
transition metal atoms in a distorted octahedral crys- 
talline electric field, the value of the Van Vleck suscep- 
tibility x VV increases as the deviation of the g value 
from the free-electron value of 2 increases. Thus from 
the g values determined by ESR for HP-(VO)2P20y in 
Table 0, one would predict that x^ V > xt W Z X^ V for 
this phase. This expectation is borne out by the values of 
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X" in Table VI for a high-quality single crystal of this 
phase. This agreement further supports our conclusion 
that the x(^) data are consistent with the presence of a 
single type of S = 1/2 AF alternating-exchange Heisen- 
berg chain in this phase. The powder average of our x VV 
values in Table || is close to the value 6 x 10~ 5 cm 3 /mol V 
estimated froni_x(T) and 31 K(T) NMR measurements by 
Kikuchi et alM for AP-(VO) 2 P 2 7 . 

Additional confirmation of the two-chain model for 
AP-(VO)2P2C>7 is the agreement we find between our 
predicted one-magnon dispersion relations in the chain 
direction for the two chains with the results of inelas- 
tic neutron scattering measurements at small and large 
wavevectors. In the intermediate wavevector regime, our 
calculated dispersion relations of the respective chains 
predict that they should cross. To our knowledge, there 



are no relevant inelastic neutron scattering data yet with 
which to test this prediction. With the caveat given in 
the next paragraph, our final estimates of the intrachain 
exchange constants and of the spin gaps of the respective 
alternating-exchange chains in the two phases are given 
in Table |[ where the error bars on each quantity take 
our mean-field modeling of interchain interactions into 
account. 

By fitting the experimental dispersion relations per- 
pendicular to-.the two chains in AP-(VO)2P2C>7 of 
Garrett et alLlL by the theoretical predictions of Uhrig 
and NormandE-3 which incorporate the influence of inter- 
chain couplings, both of the couplings J a and J x were 
found to be small but nonneghgible. In addition, the 
theoretical dispersion relationsa show that these cou- 
plings change the spin gap from that of an isolated chain 
with the same intrachain exchange constants, whereas 
our modeling of the experimental x(T) data including the 
influence of the interchain coupling in a mean-field ap- 
proximation implicitly assumed that the interchain cou- 
plings do not change the spin gap. Using our interchain 
exchange constants obtained from fitting the experimen- 
tal dispersion relations perpendicular to the chainstll us- 
ing their theory and using our intrachain exchange con- 
stants obtained from modeling the experimental x{T) 
data, the spin gap of each chain was calculated using 
their theory to be significantly smaller than the actual 
spin gap for each chain. Thus the intrachain exchange 
constants we obtain from the mean-field treatment of the 
interchain coupling should be considered to be effective 
values within this model. An improved evaluation of the 
exchange constants from xC^O data will only be possible 
using a theoretical expression for x(T) which incorpo- 
rates the effects of the interchain couplings on the two 
spin gaps. 

The spin gap of HP-(VO)2P207 obtained from analyz- 
ing x{T) data using the low-T approximations X*{t) = 
(A/Vi) exp(-A*/t) [Eq. @] is found to be different 
than obtained using the above high-accuracy theoretical 
X*(t,a) function for alternating-exchange chains to an- 
alyze the same data. For example, from a comparison 
of the spin gap obtained previously for a powdp. sam- 
ple of HP-(VO)2P207 using this approximation! 3 !] with 
the spin gap we obtained from a fit to the same data 
set using the accurate x*(^ Q; ) function, we infer that 
the error involved in determining the spin gap using this 
low-T approximation is about 50 % in this case. Similar 
discrepancies have been found previously when analyz- 
ing x(^) data for ID spin systems in a similar way. In 
the compound SrCu203, for example, the spin gap of 
the S = 1/2 Cu +2 two-leg ladders within the CU2O3 
trellis layers obtained by fitting x(^) data up to tem- 
peratures T ~ A/fce using Eq. ( |7a| ) (assuming that A is 
an independently adjustable parameter) yielded A//cb = 
420KJ13 whereas inelastic neutron scattering measure- 
ments on this compound yielded A/fee ~ 380K.C3O On 
the other hand, we found that the low-T approximation 
X*(t) = [A/t) exp(— A* /t), in which the power of t in the 
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prefactor to the exponential is modified, can yield much 
more accurate values of the spin gap. 

Our AF exchange constants in Table | along the 
alternating-exchange V chains in the two phases of 
(VO)2P2C>7 are of the same order as the nearest-neighbor 
exchange interactions estimated experimentally^ and 
theoreticallyo between the V ions in the two-leg ladder 
compound MgV2 Q». Jaut-are much smaller than the value 
of 660-670 K founcflBEHla for the V-V coupling across 
the ladder rungs in isostructural CaV205. Korotin et al. 
have inferred theoretically that the large differences be- 
tween the exchange constants in the latter two com- 
pounds arise from the stronger tilting of the VOspSauare 
pyramids in MgV20s as compared to CaV205.c3 The 
conventional empirical rules for estimating the strengths 
of nearest-neighbor superexchange interactions in oxides 
are strongly violated in CaV205 and also in cuprate spin 
ladder compounds, as extensively discussed in Ref. ||. A 
similar analysis of the exchange coupling strengths in the 
two phases of (VO)2P2C>7 would be informative and per- 
haps quite relevant to a more general evaluation of this 
issue. 
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